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ABSTRACT 


monalhe Processing for Antisubmarine Warfare is a short 
course in electrical signal processing fundamentals and their 
maieireacions in the field of antisubmarine warfare. It con- 
mans al introduction to Fourier transforms and their proper- 
ties, sampling and quantization, filters and bandwidth 
requirements, random signals and noise, and an introduction to 
four types of processing equipment; the DELTIC, energy detec- 
tOrs, correlation detectors, and beamformers. Course objectives 
meemwoiven in terms of specific questions which a person com- 
pieting the course should be able to answer. The course text 
e@deiliustrative material is contained in the appendix to the 
merce fic course 1s designed to be presented in the Fieet 
jeethe personnel involved with the operation and employment 
of detection equipment to provide them a better understanding 
Samene Operations accomplished by their equipment and to develop 
in them a better appreciation of the problems and limitations 
associated with signal detection in the antisubmarine warfare 


environment. 
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INTRODUCTION 


Signal Processing for ASW was developed as an introduc- 
mory course intended for the Naval Officer who has had the 
traditional training in ASW hardware operations and ASW ad- 
Mimtrstrative requirements. However, it is applicable to any- 
em@emtnvolyed in ASW - officer, enlisted, or civilian - who 
does not have an engineering background. 

itemce are several tields that are of particular interest 
Meme AoW Specialist: acoustics, Signal processing, environ- 
mental sciences, and probability/statistics. This course 
maresses only the signal processing aspects of ASW. Short 
meses COVED “Che otnerm Lhree areas are aiso under deveivup- 
Memt at the Naval Postgraduate School. Other theses relevant 
lgomtnese areas are KING and SWIFT, 1975; DECKER and SOLLEN- 
Miiveek, 1975; and SMITH and VICE, 1975. For anyone who takes 
all four of these courses, it is recommended that the signal 
processing be taken last, because the information contained 
in the other three will enhance the student's ability to 
obtain an understanding of the capabilities and limitations 


Ssicnal processing. 





COURSE OBJECTIVES 


mite wOwcrd bl eobjJective of this course is to familiarize 
Mme Student with the principles of signal processing used in 
mow sensor systems. The student who successfully completes 
this course will be knowledgable in the specific areas which 


moeerow to the degree indicated. 


Peo LGNAL PROCESSING THEORY 
Understand the theory of signal processing sufficiently 
melt to do the following. 
Meese ene relationship of the Fourier Transform to 
meee time domain" and the “frequency domain". 
MEEGlwen ede rOurLer Transform pair, identify the opera- 
mons they express. 
3. Given a square wave and a rectangular wave: 
a. sketch the correlation of one with the other; 
PS coLreimtEMe CONVOlULION Of One with the other; 
c. transform the given waves from the time domain 
to the frequency domain, multiply them, and apply the inverse 
franstoOrm to return to the time domain; 


dweeGomedre tne results obtained in a, b, and c. 


Bee OLGNAL PROCESSING FUNDAMENTALS 

Be ee teemer tt novledeab le Hest eprocessing funda- 
temeals tO do the following. 

1. Show the difference between a digital and an analog 


Saenal . 





Eee Urscuss Now morse limits quantization in the sampling 
and quantizing of a signal. 

fPeeerelate oampling Theorem, Nyquist Rate, and “aliasing". 

4. Given a vlot of signal and noise as power density 
menrsi@s trequency, show how filtering increases the ratio of 
Signal power to noise power. 

MeCiveiinde plotwot filter response versus frequency, 
identify a high pass, a low pass, and a bandpass filter. 

6. Relate filter bandpass, integration time, and 
mecauency resolution. 

fee Uexine the statistical characteristics of a random 


Signal and relate them to the signal voltage (or current) 


Pucumiheedssumptions mecesSary to apply statistical 
fmeenmiques to the processing of random signals. 

[ecewrnast the cross Correlation of correlated signals 
feenecne Cross correlation of uncorrelated signals. Show how 
the correlated signals are enhanced. 

Hex olaimethe three effects of filtering gaussian 
white noise. Discuss the relationship of these effects when 


Processing Signal and noise. 


C. DELTIC (DElay TIme Compressor) 
Bemdnbcem to describe the DELTIC, including: 
Po new the DELEIC functions; 
2. the benefits of using a DELTIC in the processing of a 


Siena |: 





Me TescimencespCeltications of a particular DELTIC, cal- 
Embate the integration time required and the frequency reso- 


meevon of the output. 


eee ENERGY DETECTION 
Pesecurpye energy detection by explaining: 
1. how the system determines that a target is present; 
Peeuiewe the threshold setting affects the probability of 
femection, P{d), and the probability of false alarm, P(fa); 
emetic wSi@niticance of the direct relationship of a 
P(fa) with a given P(d) as represented by a Receiver Opera- 


ting Characteristics (ROC) curve. 


Pee CORRELATION DETECTION 

be able to: 

ieeeccccripe the method used to discriminate between the 
returning echo and noise; 

CSG ume) crtect Of doppler on correlation; 

3. explain how doppler is determined; 

ecco oemtpe How pimse Length efrects doppler resolution; 

eedeserime Mow pulse length effects range resolution; 

° 6. list the advantages of using an FM slide; 

7. state the advantages and disadvantages of using 
Pseudo-Random Noise pulses for search and tracking; 

8. given a two-hydrophone array being used with a corre- 
Baumtoen detector, sketch and explain how the signal direction 
is obtained, why a direction ambiguity exists, and how the 


ambiguity may be resolved; 





Pe expldaieune lima tations on the type of signals which 
fieweoe processed passively by a correlation detector, and 


why these limitations exist. 


F. BEAM FORMING 

mee Describe the method used by the beam former to obtain 
emrection. 

Z. Name the type detector used with a passive beam 
former. 

53. Given a hydrophone array in a beam former configura- 
meen sketch the output of the beam former for a signal in 
the beam and a signal out of the beam. Show by comparing the 
Outputs the importance of setting the threshold to determine 


Mmeareet present in the beam. 
ii ou DENT 


meee ted that the Student will be currently in- 
fomved in ASW. Prior courses in probability, college algebra 
emeateiivs, and basic electricity are helpful but not re- 


quired, since the course is self-contained. 


THE INSTRUCTOR 


There are three basic qualification needed by the instru- 
HemetOr this course: a background in electrical engineering 
that includes the fundamentals of signal processing; a know- 
Peocmor CUrrent practices in ASW; the ability to,teach. The 


Naval Postgraduate School has faculty available with these 





feeiLications, and therefore should be considered as a logi- 
Soe Ource. NPS graduates may also provide a source of people 


with the requisite qualifications. 


COURSE DEVELOPMENT AND PRESENTATION 


This course was designed to provide a thorough understand- 
niet the principles involved in the application of spectrum 
Bey Sis to ASW signal processing, and the operation of the 
Mme. energy detectors, correlation detectors, and beam 
formers. These processing methods comprise the core around 
fimeeme the U.S.NAVY's ASW sensor systems are built. In pro- 
Seeding toward the above goal, the course follows a logical 
Saeemoach Of building from the basics of Fourier Transforms 
ememrourter Transform Properties. Followins those two sec- 
imomis 1s a2 section on Sampling and Quantization and a section 
mierrlters and Bandwidth Requirements. Thesé latter two 
Peervons may be interchanged, but both are required for the 
weerl0n On Random Signal, Power Spectral Density, Baie Noise. 
fimesmethods of processing are presented last; in order Ene 
PiemeDELTIC, Energy Detection, Correlation Detection, and 
Beam Forming. Where appropriate, specific hardware has been 
mentioned, but it is emphasized that the course is primarily 
devoted to principles and methods that may be applied to var- 
ious systems employed on either air, surface, or subsurface 
fea trorms . | 

Biticmeolmoceismcdestanead tO be taught in the Fleet. 


Although it requires no outside references, each student 





should be provided with a complete copy of the course materi- 
al for his (her) use and retention. The course objectives 
mmemevery Specitic. They are designed for use by the instruc- 
tor and student as a guide for measuring achievement: they 
Mmicmeate tO tne Student what he should learn; they should 
serve the instructor as the basis for examination questions. 
Mmimonliy training aid required by the instructor is a chalk 
board, as students may follow along on the illustrations in 
Mmmemcourse material. It will be helpful if 35 mm slides of 
eeeeetene illustrations are available. Overhead projectors, 
although common at shore stations, are not common on board 
meepeewhnereas a 55 mm slide projector is normally on board 
momeotner training requirements. References to books which 
cover the subjects in greater detaii, and which have been or 
ieeamstance to the authors in the preparation of this course 
aees included in the bibliography. 

If the course is presented independently, no more than 
two hours per day should be devoted to in-class presentation. 
The student should plan to spend a minimum of one hour of re- 
view and study for each in-class hour. The student must have 
time for the information to become impressed in his memory 
Beemepceone familiar to him. If all four courses are present- 
ed concurrently (feasible in a training environment due to 
the availavility of the students time and approximately the 
same time required for each course), no more than one in- 
class hour per day should be spent on any one course. Regard- 


less of the method used, the instructor should be available 


10 





mmpoitner than scheduled-class time to clarify points on 
which the student may have questions. 
The course is designed for twenty in-class hours. The 


gn-class time requirement schedule is as follows. 


Se LON FOURS: 
imereduction to Fourier Transforms- - - - - - - - - 4 
Fourier Transform Properties- - - - - -----+-+- 4 
Sampling and Quantization- - - - - - - - - --- -- J] 
Random Signal, Power Spectral Density, and Noise- - 4 
i IOI ee Ge) nS SS Se Spe ee 7] 
bmlereay Detection- - =- - - =- - - = --+-+e*+ + -*+ +--+ - 2 
Werrelation Detection- - - - - - ---s+*+7*e-+%e7e- | 2 
Beam Forming- - - - - - - - 7 - 7 = - 27 = = = = = = J 
iiemifcmerla:r COVerca is CGUuivalcnt to that covered in 


a two-credit hour, one quarter course at NPS. The sections 


Whaech required the most time to develop were the two 
eawecvens concerning Fourier transforms and the section 

Si Random signal, power spectral density, and noise. 

imechese sections the emphasis is on development of theory 
miameoncept which proved to be the most difficult to 

feetain clearly and in a manner to preclude misunderstanding. 
The other areas were eaSier to address as they lend them- 


Peweves mucl) more readily to a descriptive explanation. 


11 





CONCLUSIONS 


imecomeeurse erfectively covers the major areas of signal 
processing applicable to ASW. Its development has been di- 
rected toward the understanding of theories and principles 
used in current ASW signal processing methods. These theo- 
ries and principles will remain, although in future proces- 
Sing methods the relative importance of their individual 
moemications may be altered. For example, the current trend 
in ASW signal processing is toward more utilization of digi- 
Memeeprocessing methods. This will increase the importance 
of thoroughly understanding sampling and quantization. 

An area which remains to be addressed is the use of 
digital computers in signal processing methods. This should 
mielude an introduction to boolean algebra, elementary logic 
and set theory, the use of truth tables and Karnaugh maps, 
be ry arithmetic and electronic logic circuits. It is 
recommended that such a course be developed as part of the 


continuing education program at the Naval Postgraduate School. 


ee 





APPENDIX 
iN CRODUGCTEON TO FOURIER TRANSFORMS 


An introduction to signal processing requires, for 
understanding of many processing schemes, a working fam- 
iliarity with Fourier transforms. To gain an understanding 
of what Fourier transforms are, what they can do, and 
how to use them, let us first examine a common type of 
Signal. The electrical power used in the United States 1s 
60 Hertz (cycles per second) alternating current. That is, 
the current changes its direction of Elow in such a manner 


as to complete a cycle sixty times per second. If we plot 


the amplitude of the current in a given direction as a 


miimetion Of time, the result is a sinusoid. 


— CS i 


AMPLITUD & 





Since this is a description of the signal as a function 
of time, it is said to be the "time domain“ description. 
imemwe plot the signal as a function of frequency, we get 


the "frequency domain" description. 
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AMPLITUDE 


en 
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pimee the Signal consists of a single frequency, there 
is but one frequency represented on the f-domain plot. 
Meteus mow Look at a different signal, say an alternat- 


ing signal of 130 Hz. The time domain and frequency domain 


pots are as follows. 
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Note the way in which the changes in the signal are evidenced, 


both in the time domain and the frequency domain. 





Hz and 420 Hz. 


Likewise, two other signals, 300 
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SOmtar We have looked at signals of a single frequency. 
We are going to look at signals made up of many frequency 
components, but first let us review the concept of phase. 
PHASE: The phase of a sinusoidal signal with respect to a 
Mmecerence Signal is the relationship between corresponding 
Barts of their cycles in the time domain. For example, take 


a signal for reference. 


\ 





AMPMITUDE 





_ 
J 


Comparing a signal to this reference we can determine 


ry 


Bree pnase of the new signal with respect to the reference. 
If the new signal crosses the t axis at the same time and 
miethe same direction as the reference, it is said to be 
Sin phase’ with the reference, or to have a phase angle of 
O°. If the new signal crosses the t axis on the down-swing 
mien ehe reference signal crosses on the up-swing, as shown 
meen lt 1S Said to be "180 out of phase'’. The reason for 
using "180° will be shown later in this development. For 
miesmoment, nNmowever, just accept the fact that a signal 1s 
folds to go through a phase angle of 360° each complete cycle. 
Riempiacceeoreaestonal can either lead or lag the reference. 


That is, the signal may cross the t-axis on the down-swing 
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before Or atter the reference does. The first case would 


im—ereading and the second, lagging with respect to the ref- 


Peiecw ns bccause Gf this, a signal iS not generally described 
as having a phase angle greater than 180° since a 180° lead- 


ing signal looks just the same as a 180° lagging signal. For 


example, a signal which is leading the reference by 270° 


looks the same as one lagging the reference by 90°, and is 


Movialiy described as the latter. 





Signals with Multiple Frequency Components: Returning to 


the signal with more than one frequency component, let us 


[oeheat the Way in which this would be represented 


iL? 





me the time and frequency domains. Take, for example, the 


mone frequencies we looked at eariier: 60, 180, 300, and 


420 Hz. If we had a signal which was composed of equal am- 
Mmeeende Components of each of these frequencies, the fre- 


quency domain plot would look like this. 


1 Omer ae 
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Semertet us look at the time domain plot of this signal, 


Starting with just the 60 Hz frequency present and then as 


memadda to the sacng! judiciously chosen amplitude and mnhase 
~--= — -— @. ~ — we we we ee ee lee -_—- _— — as -——-A ee eet 4 ot ome Oo] ] EF Ot ee 2 & Oe ak & — — 


Bemiementseot the other frequencies. First, take the 60 Hz. 


mortion. 


AMPLITUDE (A) 
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New add one-third the amplitude of the 180 Hz signal but 


eaeet it 180° out of phase with the 60 Hz Signal. 


A-> 


1&0 
Go Cs 





i ! roe 
Note the fact that the 180 Hz component was added 180° out 
of phase shows up on the frequency domain plot as a negative 
teeerivdae. ihe reason for this will become apparent shortly. 
meeenow, nowever, let us continue with the other frequencies. 


memtne signal we now add one-fifth the amplitude of ele SU 


hieeconmponent, in phase with the 60 Hz signal. 


_ 





ee es ee ee ee ce ee ee om ores ee er ee 
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ieweaad one-Seventh the amplitude of the 420 Hz componeni, 


but 180° out of phase etrometme OO Hz Signal - 
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imeewe Continue to add carefully selected frequency components 
to the signal, the time domain signal becomes more and more 
like a square wave. With 20 components, the wave looks like 


mis. 


+> 


ARAAAAAS A 
} 
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From this, one can see that by careful choice of component 
imeequencie€s, and their respective amplitudes and phases, 

that a square wave can be generated from purely sinusoidal 
semponents. If an infinite number of components were used, 
the resulting waveform would be a perfect square wave. This 
Semmecept 1S the basis on which a large part of our signal pro- 
Mmecoine theory rests. There is not time in this course to 
demonstrate it, but it is a basic theorem of Fourier analy- 


Sis that any shape waveform can be created from combinations 


Bee pure Sinusoidal signals. 


wempeex Quantities: In order to simplify the mathematics 
mivolved later in the course, we must introduce the concept 

om complex quantities. In solving the mathematical equations 
describing waves, one encounters solutions which contain terms 
foeetaipiied by v-1. As one can readily see, there 1s no real 


number which, when multiplied by itself, will give a number 


rags 





morech 1s megative. Thus, these quantitics are known as 
mmracinary:. they are, in actuality, artifacts of the math- 
Mmumetes ana are Not directly related to tangible quantities, 
however, they prove to be of use in describing the behavior 
of wave phenomena. 

meee ploOteal |) “real quantities on one axis, all "“imag- 
inary'’ quantities on another, and time on a third axis, we 


Mee a coordinate system like so: 


vores ween 


Reav 






Lrwacimany 


Ome OL the results of the solutions described above is that 
a Sinusoidal signal in the real plane (the plane formed by 
the real and time axes ., 2 Mages ne Reauat Seene- 

mam De cenerated by a 
mector rotating with a 
uniform angular veloc- 
Hty in the complex 

Tem Ad RY 
plane (formed by the 
real and imaginary 


pees). Looking at 


mis vector down the 





time axis shows: 


IY 





As this vector rotates it is also moving down the time 
feos, sO the path swept out by the tip of the vector de- 


Memepes a helical path around the time axis. 


Rear 
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Nite pment ee snould be noted that there ws another 
Petatire OF anctilar rotation known as the "radian". There 
me ct tadians in SG0°. This measure is used because it 


Piplittes the mathematics involved in solving wave equations. 


ey 





LcLrumitieeno tie Vector generator, we see that the number 
of times-it makes a complete 360° (21 radians) sweep in one 
second of time is the frequency of the signal. The rotation- 
al velocity of the vector, w, is the number of radians swept 
mt, Oy the vector per second. Thus w = 27£ since it sweeps 
Out 2m radians each 360° rotation and makes f rotations per 
Second. Thus w is known as the “radian frequency" of the 
Signal. 

itewe look at the projections of the vector on the real 
aga imaginary planes as it rotates, we find that the projec- 
muon On the real plane 1S a cosine wave and the projection 
Memecive imaginary plane 1s a sine wave. Also, we see that 
me vector rotates 360° during one cycle of the cosine wave 
on the real plane. This is where the use of 360° of phase 


pemole arose. 


Rear 






Lewewaay : 
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The vector has components in the real and imaginary 
Peanes. A mathematician, by the name of Euler, has shown 


the relationship: 





where j is taken to mean “sn the imaginary plane". 

By VoCrOnEdGartTOn. tice vector penerating our signal is 
the vector sum of the real component, cos wt, and the imag- 
miary component, j Sin wt. 

Vectrohe=TeOssat tj Sin wt 
Memsce that this 1s of the form of Euler's identity if 0 = wt. 
iavs, the mathematical description of the vector is god Be 
iis vector is known as a "phasor". 

Wemm@ave scen that GoS WL is the projection of the vector 


frp path on the real plane. Solving Euler's identity for 


mes wt, we get: 
Goceut  =-ekeal Part of [e7J wt, 


Relies! wae 


rr 5 8 ~ ° Leap ~ ° - ° 
jihe SeMVoiPeT ei Or USNs ©€ ~ Per COs Wt LS used 17 


Pimbecaluse it suggests several properties of the cosine 


mea pparent in the written notation "cosine wt". ] 


mmemaacs this mean? The + sign indicates that there are 
two solutions. This corresponds to two phasors rotating in 


mime ta and -w directions, respectively. 


REAL 
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These phasors unwind in opposite directions along the 


mtx 1S. 


: Rea 





Rear 






— ee i i ie - Lene avy 


meine at the projection of these phasors on the real plane, 


Peesee Cosut. 


Rea 


Geran 
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Bolving Buler's identity for Cos wt, we obtain 


BOE 4 i) wt 


Cos wt = 7 


Weecan see this, as the two contributions of. the phasors 
always add in the same direction on the real plane, there- 
fore, the + sign in the formula. Looking at the projection 


of these phasors on the imaginary plane, we see sin wt. 


Rear 





& — he 


Solving Euler's identity for Sin wt, we obtain 
= 3) Ne 

Sil ag = 5 g 

Meecause the two projections are always in opposite direc- 
imeons On the imaginary plane, we get the — Sign in the 
mormula. 

We have shown how a single frequency plots as a single 
mime in the frequency domain. This relates to a line for 
each phasor. But we have seen that the mathematical solu- 
tions establish two phasors for each frequency, one rotat- 
mig atc +0 and one rotating at -w™. Thus, in actuality, we 


mave been showing only one side of the spectrum, the posi- 


Mmereerreaquency side. The full spectrum of cos wt is then: 
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fame lines are both on the same side of the ™ axis because of 
meet Sign in the formula 
Cos wt 


Mire full spectrum of sin wt is, 


ure Se tee eee | 


er | 


Note that the lines are on opposite sides of the w axis due 
mmethe — sign in the formula 


Wt eae 
= 
By . 


j 
Sate S 


Meee also that the lines correspond to projections on the 
real plane for cos wt and the imaginary plane for sin ot. 


We have previously shown that arbitrary waveforms can 


eg 





Pesinade from combinations of sines and cosines of the proper 
miolatuaes, phases, and frequencies. How do we determine 
which frequencies and what amplitudes and phases to use? In 
feeect, we correlate the waveform with various Piccei@rc a 
We compare the waveform to a given frequency and see how well 
Beepwemacch. One measure of this is the point by point pro- 


duct of the two averaged over a period. ThisS measure, known 


as Co (nth Poulan BeOchilelent) 1s the following; 
1 
2 —s t 
é = 2 x(t) e JT Ge 
n al 
O 
aoe 


where x(t) is the waveform of interest, eJ¥’n® is the cosine 


Dave heine compared, and the integral serves to sum the pro- 
@uct over the period, Lee of the waveform. This determines how 
Meese ly the waveform matches a given frequency. In order to 
mena the components of a given waveform, we must compare the 
meeerorm to all frequencies. The "spectrum" of the waveform 


moethen the summation of C's Mien hcamPiyemenlentrequency 


jhe : 
component (ec) Nn") over all frequencies. 


-X(£) Vic, aint? 


Using this technique, to analyze a square wave train, we get 


the spectrum as shown on the following page. 








Wesfind that all the Ge equal zero except for a few. These 
marn Out to be harmonics of the lowest component. This is 
Meeause the signal is compared to sines and cosines. If the 
meomd! 1S compared to other forms of signals, the Cos KOUIEIL 
merenecessarily be harmonics. This technique of analysis of 
MPPMeteior 25 ce mc!Opca by Fouricr amd thus has become knomn 
memrourier Analysis. The fact that some of the cis are 
Mem@ative indicates that.the frequency corresponding to that 
C. is to be 180° out of phase with the nose Sve Sign C_ 11g - 
mmencies. 

Let us look at some square wave pulse trains and analyze 


meecir Spectra. 


le~— 'T, ~~! 


ea gh a 


If we solve for the C_'s and plot we get the curve as shown 


Pee the following page. 
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Miemcne period of the pulse train is increased but the length 
Seeeeach individual pulse is kept constant, the following 


Surves are obtained. 


ieee that the point where the envelope of the oes ELOSsSes 


zero remains the same, but the number of C's eine a+ 6S. 


That is, there are more frequency components present in this 


Emenal than there were in the. first. Let us now increase the 


ye 





Mmemrodnwevyen more, still keeping the pulse length constant. 


bo Y, 





—_— <— 


Mere miiat tie zero Crossings of the envelope remain con- 


meant, but the number of C's increases again. 


If we were 


Meme >crease the period to the point where there was 


mmeare pulse left, that is, T > @ 


On, soe 


, we would find that the 


number of c's VOU@eecmiitInives bUE the Zero Crossings Of 


mee envelope would remain the same. 


ket, -o4 


—s. | 
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From these examples we can see that the zero crossings 
mierval of the envelope is related to the inverse of the 
pulse width, and the interval between C's is related to 
the inverse of the period of the signal. This, aunval “2 Weve 
Meeoperties of the Fourier relations, will be discussed in 
Peeater detail later in the course. 

The form of the envelope of the C's for the square wave 
memears SO frequently that it has been given a name of its 
Same the form is (Sin x)/x, which is called "sinc" (pro- 
momnceed 'sin-cee"). 

Fourier Transforms: We have shown that as the period 
meamptestonal becomes infinite, its number of spectral compo- 
nents also becomes infinite. In this case, the spectrum of 
meeenal can be characterized by the formula of the envelope 
of the C's. Techniques have been developed for finding the 
formula of the envelope without solving for the c's and one 
Meumnod of doing this is known as the Fourier Transform, de- 


med thusly: 


oO 


X(£) -{ x(t) e J Yt at 


= co 
foeee chat this 15 quite Similar to the formula usSed for find- 
ing the C'S. Instead of giving us the individual C's how- 
ever, this gives us the formula for the envelope of the Co's 
mie the f-domain. 

If we compute the transform of a single Square wave using 
this formula, we get the following result for the form of the 


Felution: 


3, 





mf) = + Sinc x 


This is the same result we get if we compute the c's Le 
a square wave of T + ~ and solve for the formula of the 
emvelope. 

Bence we scan Lund the spectrum, or f-domain characteris- 
ies Of a Signal, given its time-domain formula, it is only 
reasonable that we should be able to do the reverse opera- 
tion. This turns out to be the case, and the operation is 


accomplished by using the "inverse transform". 


x(t) -| ce) el Oar 


t co 
Peeoubstitutaon in the £-domain formula of a signal and 
Beemeintescratinag, one can get the time domain characteris- 
tics of the signal. The two transform formulas, forward and 
Pero wedre called a “transform pair’. One finds the f-do- 
main characteristics given the t-domain formula, and the 
emer finds the t-domain characteristics given the f-domain 
mermula. Later we will see that these are powerful tools 
Mseed in Sipnal processing. 

For the moment, let us return to the cosine wave which 
memlo@ked at Carlier and see what happens if we introduce 
modulation. That 1S, Superimpose it on a carrier frequency 
cosine wave similar to AM radio transmissions. 

A modulated cosine wave looks as picturadon the follow- 


ime page in the time domain: 


SE 








where the envelope varies as the cosine wave we are looking 
at and the individual peaks are at the carrier frequency. 
What does this Signal look like in the f-domain? We know 


fee the £-domain plot of our cosine wave looks like. 





Mee carrier iS just a cosine wave of a different frequency 


meats plot looks like this. 
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Mete that both f-domain plots were centered about f = QO. 
fem, if we put the two together, i.e., modulate the carrier 


with our cosine wave, we get the following: 





ae 


+f.) -(4-f) “ £ an 


Tenet naeethie resulting Signal 1s of the form Cosine 
i mei). the £-domain plot is the cosine wave plot cen- 
meeead about the carrier frequency values instead of about 
f= 0. Thus, modulation in the time domain is equivalent 
Mmemaetranslation of the spectrum in the f-domain from cen- 
MmeececdeaboOut the f= 0 axis to centered about the carrier fre- 


quency. Let us look at a modulated ‘square wave pulse train. 


| 


owe 
eee ees oe SO 


ia ) 
Tmieowmel@em tic L-comain plot Of a square wave iS a 
sinc function. We have just seen that modulation of a sig- 
ime cranslates the f-domain plot of the signal and centers 
it about the carrier frequency. Thus the f-domain plot of 


a modulated square pulse train is as shown on the following 
page. 


3f 





where the envelope takes the shape of the square wave spec- 
mem Out 1S Centered about the carrier frequency instead of 


about f = 0. 
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II. FOURIER TRANSFORM PROPERTIES 
meme properties of Fourier Transforms are stated most suc- 
mmecty in the accompanying Fourier Transform Theorems. Know- 
Mmemee Of the theorems 1S important in the interpretation of 
Peeetra as the theorems express the relationship between 
time-domain and frequency-domain operations. 

Bice mt roduction to Fourier Transforms it was shown that 
Mmaemrrequency spectrum, X(f), of a signal, x(t), was the Four- 
mammionansitorm Of the signal. That is, X(f) = F[x(t)] and in- 
versely, rach athe ect) 1. thesesnre lations, are calved “transtorm 
mars” and are more conveniently denoted by x(t) «.X(f). The 
memmder Iransform Theorems describe the properties of various 
Fourier Transform nairs, 

Pee NEARITY (OR SUPERPOSITION) THEOREM 
mia Sienal, x,(t), See lete Oe ked a Diyeec constant, a, 
fas would happen if an amplifier with gain=a were 
Mdccdehn sdeclrcuit) and a signal, x5(t), Hs ge LEE 
(Pittcdtmameconschant.. bp. them their respective line spec- 


Pawiiealso be multiplied by the constant a or b. 


x, (t) => ax, ( t) X(t) bx. (t) 


ax, (t) + bx,(t) @- aX, (f) + bX, (£) 


Oy 





mie ar comDbDinations in the time-domain become linear combina- 


mieme, in the frequency-domain. 
B. TIME DELAY THEOREM 


Given the transform pair x(t) q@>X(f) if 
piceobena x(t), 15 delayed by Ee SSE: 
onds, as happens when the signal passes 
PanOUgiecmerpdcltIve OF Inductive element, 
to become a new signal, x(t-to), then the 
Spectrum ws modified by a frequency de- 


pendent phase shift to become Rie Jones 


———s 


1wt 


X(t-to)€E>X(f)e J°%o 


xt) |e. 





jhe translation of a signal in time changes the phase of the 
iecererum but does not alter the complex amplitude. 

Momebbuistrate the Linearity and Time Delay ITheorems, con- 
oder ad Signal which 1s a linear combination of delayed sig- 
nals: 


Age eee) 4x (t- St.) 
Mes Spectrum 1S easily written as: 


VP Ry a Die TOM a Tinea ore chute 


1.0 





G> SCALE CHANGE THEOREM 

rt has DeCiSiengiethatwa translation of the time origin 
may be accomplished using the Time Delay Theorem. The time 
axis may also be expanded, compressed, or reversed by an opera- 
miomeknown as ‘scale change"’. If a signal, x(t), becomes a 


new signal, x(at), then: 


memeex(at) 1S compressed if a 1S positive, greater than one; 


mC) = a .. > 1 x(at)= = 
= to t=¢ 


memeex(at) 1s expanded 1f a is positive, less than one; 


amet sis reversed in time 1f a is negative. 





In the Fourier Transform pair, 
] ie 
Xx (at) Sa ae 5 
la] means the absolute value of a, that is, the numerical value 
Only is used, the sign is dropped. 


The scale change theorem expresses the property of recip- 


rocal spreading. If the signal is compressed in time by the 
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PaeeoY 2, its spectrum is expanded in frequency by 1l/a. If 
ieesignal 1s reversed in time (a = -1l) its Spectrum is 


) = X(-f), which indicates a reversal in phase. 


X(t) aes X(é) 
t=0 ad > 


_, exe) 





The expansion or compression in the time axis of a signal 


occurs most often in the playback of recorded signals. 


D. FREQUENCY TRANSLATION (MODULATION) THEOREM 
| ; 


Piiesipmal. xlt), 1S muitiplied by eJ%c , 
Loma Peeunlmn., A(t), will be translated in 


eC Ncvweny +t owe Mathematically stated: 


x(tjet¥ct esx (£-£,), 


X (t) 


We 


xe?’ 





Le 





iemis more common to multiply a time signal, x(t), by the real 
part me el%c’ ice yecoss-—. Prom previous work with 
miber’S theorem, it 1s recalled that 

x(t) cos wet = 1/2x(t)(eI@c® + e J¥c*y, 
Taking this result and applying it to the Frequency Translation 


iieorem yields the Frequency Modulation Theorem: 


aC Oita ly er =) + 1/2X(£+£.) 





E+Fe 


See Signal x{t) is multiplied by cos wt, its spectrum 
Meets Shifted up and down in frequency by an amount f£,. This 


a _ 


nT ~~ = -—- 
i fdadtil Wa 


Meee Same resuit noticed eariier when an RF pulse ¢t 
meeiea Dy multiplying a rectangular pulse train by cos wt. 
mmenerea then, and confirmed by the modulation theorem, time 
@emwain multiplication becomes translation in the frequency 


domain. As an example, .a rectangular pulse of amplitude A, 


period t, centered at t = 0 and duration Tt is x(t) = AN(t/T). 





miemwspectrum 1S recognized as a SINC function, or 


MCE) = SA oUNe £r . 
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mOW, multiply tnat single pulse by cos w,t: 


Metye= All(ty/1) Cos Wet. 





By the modulation theorem, the spectrum is: 


Ce) ZesING (f-£ojir + 1/2 A SINC (£+f_)c. 





This is the resuit that is used extensively in single-sideband 
and double-sideband HF communications. 


E. TIME DERIVATIVE THEOREM 










THE tha Ceorivatixc 
fy Lasaw Mewes ah AeA v 9 


clgiak x(t), the result is multipli- 


cation in the spectrum Dy jw: 


me PGE) Cae WCCO 









Multiplication 
Pieee He dtue ley 








fF -—» 


Lydy 





Mieaetore, differentiation enhances the high frequency compo- 


emis in the spectrum. 


feet ME INTEGRATION THEOREM 


a 
Taking the integral [ Oe 21 Sipe eryh x(t!) 


- € 


Pome henteco miktiplying by = (di- 


viding by jw) in the frequency domain: 


t ' i i 
[xc ) INE ee a (Cre Ie 
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Integration, then, suppresses the high frequency components 








—_ —- mae 


but the low frequency components are unaffected. 
fee Convolution 

Take two signals in the time domain t', for example, x(t’) 
and y(t’). love elapiethe y(t } cia romeo Otniereside of 


the axis and it becomes y(-t') as illustrated below. 





' ( 
{ a ae Be | 
x(t) : — led rer F x(t) 
Fao ay 
— Eyes: + r a oe 
<—. t —> | 1a-t{ ———> ’ 
; 120 est al tz0 


t--t 
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In order to make a comparison between the two signals, multi- 
ply them together at every point in time t’ from -© to +o, 


Mathematically, this can be expressed as 


2 ef x(t) y(-t')dt’ 


meme above case, anywhere in the t' domain where et) 
amots, y(-t ) = 0, and where y(-t') exists, Ge) = 0. There- 
Meme. their product is zero. 

Since the signals are separated by time t, if a time delay 
of t is introduced making y(-t') become aera the signals 
mime overlap. If t, the time delay, is varied from -~ to +o, 
the meet ) Soiawewrll slip” across the signal x(t'). The 
Meewlt 1S expressed as a function of time, t: 

+00 


ZiGe)) -f mit yoy(tet jdt’ . 


co 
The above integral is known as the "convolution integral" and 
Ht Expresses convolution mathematically, i.e., the area of 
Mresproduct of the two signals. A more convenient notation 


fer convolution is: 
z(t) = x(t) * y(t), 


A graphical interpretation of convolution, as shown on the 


following page, may help to see what happens to the signals. 
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Two signal "flipped" 








5 
i 
{ 
j 


Gt? au 1 xe 
¢ t 
1 ’ 
-——_ cee —_ =~ a ga =e =_-- 


—t —> | 
Vere est 





cee 
fiessiaded area on the | | Z(t) 
feiemp.1ot Of both signals a ; 
feeeesents the area of their ) | —_ 
meoenmon product" for a partic- : 

meteevalue of time delay, t. Ra _A 

The snaded area on the z(t) _ ie | 

Peete represents the summation 


Smee area of their ‘'common 





mmeoauct'’ as the "flipped" 


feeemial, y(-t), is "slipped" 
| 


KW 


Meross x(t). 





ae CO 6G) 
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Micmewmniiiitonegm tt. Cc. tne “flip and ''slip', of two rec- 


tangular pulses produces a trapezoidal pulse with a base of 


length ia Oo. Convolution of two identical recrancular 

meses produces a triangular pulse with base = 2t and height = 
” 

fT. 


' ‘ ! ay ates 





| xt) xt’) 
A — . 
JRE 2(t) = XO * YW ae: 
Tt: Hitec nS tance , | mi | 
® bee -__> 
a ead Pere eet = x(t) 


Now that convolution in the time domain has been explained, 
Memedoes it relate to the frequency domain? There are two con- 


Volution theorems which relate the two domains. 


| ae) a : Y(f£ | 


x y(t) & X(L) ) 


Convolution in the time domain trans- 


Peses eOmmuUlcC@plteation in the fre- 


quency domain. 





X(t) - y(t)@>X(f) * Y(£) 


Miter iicatnon im the time domain 
PEM nOscomGOmecolyvOliLion in the fre- 


quency domain. 


For example, look at the now familiar square pulse shown 


on the following page. 
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H. Correlation 

In statistics, when it is desired to know how closely 
Simienwdaistribution resembles another, the two are overlayed and 
mmemareas they have in common are multiplied. This "common 
memeuet' represents the degree of similarity, or how well one 
Memerioution correlates to the other. In signal processing it 
memorten necessary to compare one signal with the same signal 
moameen Occurs later in time or with another Signal. Signals 
fey) in time in contrast to distributions which are constant. 
Mmmesire. iherefore, in order to completely overlay one sig- 
nal on another, displace the signal in time by tT units, then 
tt to “slip” the displaced signal across the one with 
Mimech it is desired to correlate. With each different Tt used, 
mie common product” will be a different value. In order to 
obtain one number to represent the value of correlation, sum 
mime common product"' over all values of delay, t, and aver- 


age it by dividing by time T (or multiplying by ae This 
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Seeeiiccethe. €Orrelation function" R(t) and is described math- 


eMmatically as: 


R(t) = = e a Klett dt. 


This is known as a "time average" and is denoted by 


R(t) = <x (t) (am) 


eee bans 
where the < » (brackets) indicate Too F prow MCU dt. 
meee signal is periodic in T, it need only be averaged over 
Seeperi0d, T,, instead of all time. The correlation function 


then becomes: 
i. 
R(t) = 1 f Mitevextter)at « 
To To 


iaeeeas more precisely referred to as the autocorrelation func- 
Mmmm ot 1S, the Correlation of one signal with itselt. It 
ieee desirable, as mentioned earlier, to correlate one sig- 
iemeeex(t), With another signal, y(t), displaced Tt units in 
minesto become y(tt+t). This results in the cross correlation 
mmc tion, Rey (7) = Gale) y(t+t)> We WilcmaeOGOTGelatlonm 1Une - 
mene is mOre commonly used. Where "correlation function" is 
mead it 1S assumed to be the “autocorrelation function". 
homececraphical representation, the familiar square pulse 
is used. This, as shown on the following page, represents the 
meutt Of the Summation of the “common product" as the delayed 
Meese 1S “Slipped” across the other pulse. This initial part 
Memeorrelation is like the previously described Sen olution 


except that the signal is not flipped. 


SO 








Spies 


The next step is to take the time average of the summation of 
Mmemeonmnon products. It 1s seen that, as in convolution, cor- 
relation increases the width of the function by es Ts; Oiler 


Mmeene case of a square pulse, from T to 2T. 


te 


. | 
| 
i 
3 
| | | Sum of 
= tie common products = = oy et 


Momcime average, take the 


area = 1/2 bh, or iL 2 (CRI on ; | 


and divide by time, 2T, 2. | 
i 


M/2(2T)b~ _ db 
— 2; COM 7 Z i eS i 


What does the correlation of a sinusoidal signal look like? 
meeers See on the following page. 





When both signals overlay 
eempletely, the value is 
boree and positive. 


A small shift will give 
some negative area but the 
mesmit 1s Still positive. 


A 90° shift makes the > 
meeatcive and negative areas 
equal and the resultant val- 
Mien zero. 


A shift of 180° results in 
Gmeeentire area being nega- 
mere . 


mez 70 the positive and 
negative areas are again 
Semtal and the resultant 
mate 1S zero. 


feeooU! ,a complete one cycle 
mrt, the signal is back in 
miase and the result is a 
maximum positive value, 








The correlation function of a sinusoidal signal is another 
Sinusoid, and, if the time average is taken, has an ampli- 
bé | 
mace of >: 
Up to now there have appeared many similarities between 
convolution and correlation. The graphical Piverpretat lon 
has also shown this similarity but much of that is due to the 


choice of symmetrical signals for illustration. Now, look at 


the side by side convolution and correlation of a saw-tooth 


Signal. 
Convolution Correlation 
eebrp then 'siip" Slee! wevailye 


a 
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m@emtinal results are indicated in the next diagram. 


When "time Averaged" 
becomes: 


So it is seen that if an asymetrical signal is used, the finai 
result obtained by convolution may be quite different from ele 


final result obtained by correlation. 
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II]i. SAMPLING AND QUANTIZATION 


Meche real world, most of the Signals of interest to an 
operator of a detection system are of a continuous, or analog, 
mere. On the other hand, most of the more sophisticated 
Meeecessors in use today are digital in operation, and thus can 
moetise analog inputs directly. Some means of digitizing the 
emaerog Signals must be used. 

momarcitize a2 sienal, one must describe the value of the 
Signal at any given time by a digital number. This requires a 
measurement and conversion of the signal value to digital units. 
The question of how accurately we must measure the value of 
miemciconal in converting it to digital form must be considered. 


( eres Nees Eaguere 


oMiowmenNo sled. LOT eExaliple: 


t-~ 


This signal has a maximum value 
mor 0. since we must convert the 


miemal to a set of discrete values, 


Oil Seanasd minimum val- 
COmeilmucus Values of the 


we must then decide how. 


many discrete values we need in the range from 0 to 3 to de- 


m@emaoe the Signal adequately. 


Ifewemdectac to use four levels, 


me 1S, digitize to the nearest integer value, the signal 





mm have the following appearance. 





me we double the number of values and digitize to the near- 


@se half integer value, the digitized signal looks like this. 


! \ 


8 
we wb Re ete a 


eae 


Note that as the number of values used ae Gee Zine eke 
meet increases, the digitized curve more and more closely 
ieenoximates the original signal. The values used in digi- 
fmezineg the signal are known as quantization levels. In the 
menor a noise-free signal, one can digitize the signal to 
mi eaecree of accuracy by increasing the number of quantiza- 


mon levels. 
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The case of a signal with noise is different, however. 
In this case, there is a random fluctuation of the signal 
value because of noise even when the value represented by 
the signal itself is unchanging. The range of these noise 
fluctuations determines how accurately we can measure the 
value represented by the signal even if we have a perfect 
measuring device. To illustrate this, consider first a sig- 


mepewithout noise and then one with noise. 
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WITH NOISE 


We observe that in the noisy Sanco Ene a ailuemon june 
Signal fluctuates about some mean, or average, value within 
some limits. If the noise is truly random in nature, the 
mean value about which it fluctuates is then the actual 
value of the signal being represented. At any given instant, 


any value of the signal between the upper and lower limits 
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@e the noise fluctuation could stand for the actual value 
being represented. The accuracy with which this signal 

can be measured is then limited to the range of fluctuation 
due to noise. For example, suppose the signal of interest 
memaeconstant 3.0 volts, but because of noise the signal 
Fluctuates between 3.5 and 2.5 volts. The most accuracy 
with which we can measure this signal is to the nearest 
integer volt. 

Once we determine how many quantization levels are 
necessary, we mjst then determine how often we need to sam- 
ple the signal in order to reconstruct it accrately. 
Obviously, if the signal never changes value, but remains 
Seeetant, only one sample is needed to reconstruct the sig- 
nal for all times. If, however, the signal is not constant, 
we face a real problem. 

How often do we have to sample a given signal? Well, 
we know that one sample is not enough if the signal is not 
constant, and we also know that if the signal is sampled 
Continuously (i.e., in an analog fashion) we can reproduce 
it exactly. To sample a signal at extremely short intervals 
1s very costly and difficult to accomplish. What we are 
looking for is a compromise, a sampling rate high enough 
that the Signal can be reproduced from the samples but not 
any higher than necessary in order to keep cost and equipment 
complexity down. 


We can see intuitively that there must be some relationship 


between the rate of change of the signal and the sampling rate 
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necessary to be able to reproduce it. If the signal changes 
wemy Slowly, we need sample it only infrequently. If it 

changes rapidly, however, we must sample it frequently enough 
that it does not have time to make several changes in the inter- 
val between samples or we lose these changes and are thus unable 
memreproduce the signal. If the sample intervals are too long, 
the signal is said to be "undersampled". This condition, in 
@adgition to causing us to lose some of the information in the 
Signal,can also cause other complications. Perhaps you have 
Geerced that in some of the old cowboy movies the stagecoach 
Mitece!S appear to turn backwards at times. This is a case of 
midersampling which illustrates one of the other effects called 


Saeiasing'’. To explain this, let us lock at a sine wave. 
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Ii we sample this at a frequency less than that necessary for 
reproduction of the signal, we can get a series of sampies which 


pive the impression of a sine wave of lower frequency. 





This is a problem which must be addressed when analyzing un- 
mann Silenals for frequency components since aliasing can give 
Selmerous results if the analyzer sampling rate is not high 
emewen, this will be discussed further in a minute. In the 
meantime, let us look at the sampling process. 

In the time domain, we can visualize the sampling process 
ideally as the multiplication of the signal by a function known 
as the "ideal sampling wave", which is identically equal to 


zero except at the sample times I[T., 2T_, etc.]. 
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| “IDEAL SAMPLING WAVE” 
The duration of these spikes, known as "delta functions", is 
miednaitesimal, approaching zero, while their amplitude approaches 
mimity, and they are defined such that their area (0 x ») is 
Samael to 1. When the signal of interest is multiplied by the 
ideal sampling wave, we get the product of their areas at each 
fot PT. Since the area of each delta function is 1, and ie 
meath is 0, the product of areas is equal to the value of the 
Seemal at that time. Thus the products of the two functions 
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het us look at the frequency domain to see more clearly 


what happens when a signal is sampled. 


Assume that the signal 
meee Sampled is strictly "bandlimited", that is, 


it has no 
meegwency Components outside of some limit. 





me ) uy he. 


Meas the bandwidth of this signal. 


Although the mathematics 
m@oerather involved, 


Me Ganenewne shown that the adeal sam- 


pling wave spectrum appears as in-the next diagram. 
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We know from the modulation theorem that modulation of a func- 
meme causes its spectrum to be shifted and centered about the 
modulation frequency and otherwise does not change it. Thus, 


the spectrum of the sampled signal appears as follows. 





meotiecnis we can see that for strictly bandlimited signals, if 
- Meee there 15 no overlapping of the spectrum components 
fmemeso the spectrum of the original signal is exactly repro- 
mieed. iif the sampled signal is then passed through a low-pass 
M_imeeet which has the property of allowing only those frequencies 


Meeeowe 1tS cutoff to pass, we can exactly recover the original 


Signal. 
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It foliows from the above that the sampling frequency must be 
Mewes 1.€.5, Sreater than, or equal to, two times the highest 
meeamency component of the signal if we are able to reconstruct 
the signal from the samples. If it is not, the Sle © e nuiienoss 


maeesampled signal has overlaps which introduce errors into 


mmemneconstruction of the signal. 





{-» 


—— 


Mitese Overlaps are the cause of "aliasing" and other incon- 
meeeemecieS in the reconstructed signal. This result is formal- 
ized in the Sampling Theorem, or theorem of uniform sampling: 
Piea stenal contains no frequency com- 
ponents for |f| > W, it is completely 
described by instantaneous sample values 
uniformly spaced in time with period 


rh = aon 


mies rate of sampling, f = 2W, is known as the Nyquist Rate. 
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fears the absolute minimum sampling rate from which a signal 
mame reconstructed under ideal conditions. 

‘Until now we have considered only ideal sampling and recon- 
meemerron. in practice, the ideal is rarely, if ever, realized. 
Practical sampling differs from ideal sampling in three obvious 
meopects: 

(1) Practical sampling waves are not composed of delta 
functions but have a finite width. 

Pee Practical filters are not ideal. 

(3) Real world Signals normally sampled are not strictly 
Mama limited. 
Memeractice these differences are of small enough significance 
that for all normally encountered cases only the last is of 
meet importance. 

Since most real signals are not strictly bandlimited, there 
Mie be SOme overlap in the sampled spectrum. If the spectral 
components outside some nominal value W are negligible, however, 
mieesienal can still be adequately described for most purposes 


byesamples spaced To <. 1/2W. 
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IV. FILTERS AND BANDWIDTH REQUIREMENTS 


ie ocean 1S an extremely noisy anil Rommene iio eo ueal 
to search for submarines. Noise from waves, wind, breaking 
surf, shipping, sonic mammals, fish and crustaceans, 
feeeomrc aCtiVity, rain, etc., iS present to some degree 
meal times. How do we go about sorting out these dif- 
mrent Sources and detecting submarines? One way is to 
use filters to eliminate as much of the noise as possible. 
mempccmeral, there are three classes of filters relative 
to spectral response: iow-pass, high-pass, and band- 
mse tilters. Ideal low-pass filters allow frequencies 
iow their cutoft frequency to nass, and screen out 
ai above cutoff. Their spectral response is shown in 


the following diagram. 
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Hiueh-pass filters allow all frequencies above their 


@otott to pass and screen out all below. 
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Pamad=-pass filters allow all frequencies between their 
mere ang Upper cutoff frequencies to pass and screen 


omt all others. 





bitters are characterized by their upper and lower 
Cutoff frequencies, or their center frequency and band- 
Maes) (Or bandwidth), B. The bandwidth is defined as the 
@assband width measured in positive frequency only, 


Phat is, the difference a - ft) = B. 
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Miecourse, no réal filters have the sharp cutoff 
EMaracteristics of the ideal filter, so it is customary 
Mmemaetine the bandwidth of a real filter as the band- 
Width between the points where the response has dropped 
Segoe itrom the maximum. These points are known as the 
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Peetonal which 1S passed through a band-pass filter 


meeears aS in the next diagram. 
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As indicated earlier, filters eliminate as much of 
mime moise as possible. How does the filter actually 
improve detection capability? To answer this question, 
we examine the way filters help improve detection in a 
large class of detectors, known as energy detectors. 

We discuss the way these operate in greater detail a 
little later, but the general principle of operation 


~ 
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memenat tne detector senses the total power (or energy) 
Meeeetic in the input, and compares this to a "threshold" 
Waaue. if the input exceeds threshold, the detector 
Mremreates the presence of a signal, and if the input does 
ieee xecced threshold, it indicates no signal present. 
mrmonracr tO see better how filtering helps to detect 
feeieis im noise, consider the “power spectral density" 
MmeoteOr an input. The power spectral density is the power 


Mmeaesignal that is carried by each frequency component. 
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ares total power in a signal then corresponds to the 


mmmerral Of the power spectral density over all frequencies. 


The power spectral density thus gives an indication of 
which frequencies carry most of the power in a signal, 
and the relative amount of power at a given frequency 
mempared to that at others. The power in a given freq- 
momey band 1S then proportional to the area under the 
mewer spectral density curve between the lower and upper 


meets Of the band. Comparison of the power present 
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in a signal to the power of the background noise present 
miieates the detectability of the signal. Consider 

a signal against a background of "white noise" (white 
Memse 1S defined as having equel amounts of all Speeenal 


memponents) to see the gains from filtering. 
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We can see that if the detector examines the entire 
frequency spectrum and senses the total power present 
to compare to a threshold value, the presence or absence 
Meee Sipnal does not affect the input significantly. 
Omethe other hand, if we filter the input and only allow 
Memarrow band of frequencies about the frequency of interest 
memene input to the detector, then the presence of a signal 


will substantially change the amount of power in the input. 
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mins, in effect, raises the signal-to-noise ratio, 
and thus allows us to detect weaker signals than we could 
Memect before filtering. In the case of non-white noise, 
Piemamprovement in Signal-to-noise (S/N) ratio can be 


even greater. 
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Ambient noise in the ocean is an example of non- 
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Energy detectors integrate the energy in the input 
Over some time period, some averaging the energy and 
Some not averaging. The value obtained in either case 


1S compared to the threshold, which is set at the value 
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corresponding to ambient noise when the signal of interest 
is absent. Therefore, one must determine how long to 
integrate the input before comparing to threshold. How 
iecnis determination accomplished? We saw in the section 
on Fourier transforms that there was a relationshiv 
between the period of a signal and the separation between 
the frequency components that make up the signal, namely, 
the separation between components is inversely related 


Momene period of the signal. 





Thus if we assume that the integration time we choose 
is the period of the input, and that the input exactly 
repeats itself every period ad infinitum, then the minimum 
resolvable frequency difference between components of the 
input is inversely related to the integration time. The 
longer the input is integrated, the more accurately the 
frequency components present are determined. If we 
desire to resolve frequencies 1 Hz apart, we must integrate 
meeeat least 1 second, for frequencies 0.1 Hz apart, at 


least 10 seconds. 
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DiICuier scams tderation in detection system design 
meee fact that most common filters in military appli- 
cations today are composed of reactive elements, such 
eo, LC/RLC circuits, crystals, or transducer elements, 
which have finite time constants associated with them. 
What does this mean? It means that these systems are 
eeeoned to resonate at a particular frequency or over 
meoeand of frequencies, and to block or short out others. 
mene reactive elements require a finite time to "ring 
Mgeeeor eStablish resonance. Thus, systems do not respond 
mi@eciately to inputs, and this must be considered in 
the system design. The response of these circuits has 


the following appearance. 
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fie time required for the response to rise to 63% 
of the maximum is known as the time constant of the system. 
tiie Narrower the resonance, that is, the narrower the 
miter bandpass, the eneee the time constant. Thus, 
Mmemmriemacsire to resolve frequencies to 0.1 Hz, we must 
allow enough integration time for the filter to respond 
to each frequency, and this will be longer than ine eine 


required to resoive frequencies to 1 Hz. 
3 





Pmotner Consideration in the design of active ranging 
Memmection systems is the requirement to accomodate doppler 
Moeene returning echoes. The input filter must have a 
bandwidth wide enough to pass echoes that are doppler 
shifted by the maximum amount anticipated during operation 
Mimmtne system, or the echoes will not be detected. If 
Better of sufficient bandwidth to ensure this does not 
faves a sufficient S/N ratio improvement, one may alter- 
natively use a bank of narrow band filters centered at 
mememous £requencies corresponding to different doppler 


Meets. this is illustrated in the diagram on the next 


page. 
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V. RANDOM SIGNALS, POWER SPECTRAL DENSITY, AND NOISE 


A. RANDOM SIGNALS 


Heretofore the signals discussed have been explicitly 
described "deterministic signals", such as a square wave, 
Sinusoidal wave, sawtooth, etc. In describing a signal ex- 
plicitly as sume function of time, it is assumed that its 
amplitude and phase are known exactly for all time; ieee 
past, present, and future. No real world signal meets these 
criteria, for then the signal would convey no new informa- 
mom, aS it could be predicted exactly for any future time. 
Therefore, consideration must be given to random signals. 
Random signals of interest to the Antisubmarine Warfare spe-~ 
Ccialist emanate, for example, from various machinery and flow 
noises. Noise in the ocean due to waves, wind, biological 
sources, seismic disturbances, etc., 1S a special case of 
random signal, and it is always present to some degree wheth- 
er or not a "target" is present. Processing equipment must 
be able to distinguish between the random signal of interest, 
unwanted signals, and noise. 

A random signal, unlike a deterministic signal, may only 
be described by its statistical characteristics. Its past 
Mmistory and its predicted future will be stated in terms of 


some average value and a given probability that the Signal 


76 





will be within certain limits at a specified time. This 
"given probability” is represented by a probability density 
mometi1on (pdf). 

What 1s a probability density function? A function is a 


descriptor, 1l.e., it usually describes something. Density is 


a measure of quantity contained within a specified space. 
meoodotlity refers to the chance, or likelihood, of the occur- 
rence Of an event. Therefore, a probability density function 
describes the likelihood of some quantity being contained 
within a specified space. 

How does the probability density function relate to a 


random signal? Assume a random signal, x(t), which has some 


randomly varying voltage amplitude over time. 
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Plot the signal with a vertical scale representing voltage 


amplitude and the horizontal scale representing time. 


7) 
v 
$ 4 
Ww 9 
*, 
a 9 
 S 
— 
C 





time — 


Suppose now that ]-volt represents any value between 0.5 and 
beeo; 2-volt represents }].5 to 2.49; 3~volt represents 2.5 to 
3.49, etc. From the random signal, a plot is now made with 
the vertical scale still representing voltage amplitude but 
with the horizontal scale representing the number of times 
the Signal is at that amplitude, or within that amplitude 
range. 

During time interval ], the signal was present only in 


the 3-volt range and that one block is shaded in on the plot. 


amplitude 





numher ch times the 
signal cs withen the 
ceomptctude OB 


During time interval 2, the signal is present in the l, 


2, and 3-volt ranges. Those additional blocks are shown in 
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the next diagram. 
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fieomg time interval 3, the signal is present in tne 1 


and 2-volt ranges. Those additional blocks are shaded in. 
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The development should start to become apparent. Look 


now at time intervals 4,5,and 6. The signal is present in 
the 3-volt range once, the 2-volt range twice, and in the ]- 


volt range three times. Shade those in. 


S 









y 4 

3 = 3 Oe 

co WON 

a 9 ENR NN on 

eS 4 OOM 
oO 





i tag 10 
the 
Leie 


to? fs + 5S % 


ck times 
We [Aaw 


reawge 


Numsg ar 
Signal ts 
amplitude 


@eeeinish the plot, look at time intervals 7,8,9 and 10. 


The signal 
volt range 
Eee Ii—volt 


Shade them 


appears in the 4-volt range four times; in the 3- 
three times; in the 2-volt range four times; in 
range four times; and the O-volt range once. 
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The plot above represents the amplitude ranges that the 


Signal was in for the entire time it was observed. 


The tocal 


number of shaded blocks represents the number of times the 


Signal is 


ious time 


within the various amplitude zones during the var- 


intervals for the entire duration of the signal: 


the 29 blocks account for ]00% of signal time. (Count them) 


Fach shaded block is regarded as a sample value. 


ber of sample values is denoted by n, then n = 


example. 


Le the nun 


29 in@ehis 


The percentage of time the signal is ina 
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particular voltage range can be approximated by dividing the 
number of signal sample values in that voltage range by the 
total number of sample values, n. 
v-vort twange: 1/29 = 3% 
l-Volte cange> 9/29 — 32)2% 
2-volt range: 8/29 = 28% 
Sa vObhesrmange: § 1/29 = 242 
4-volt range: 4/29 = |4% 
Total Jj003% 

The percentage of time the signal spends at other ranges 
can also be calculated. Suppose, for example, that in an ASW 
system there is a detector which detects only Signals with 
amplitudes of 3 volts or more (the detector has a "detection 
threshold" of 3 volts). To determine the percentage of tiine 
the above random signal is detectable, it is necessary to de- 
termine the percentage of time it remains in the range of 3 
Or more volts. The 3-volt range covers from 2.5 to 3.49 
volts. Take one-half the time the signal is in the 3-volt 
range and add it to the entire time the signal is in the 4- 
volt range to approximate the total time the signal is at or 
above 3 volts. 

By 247) + 4 = 7.5 Divide by n to get percentage. 
7.5/29 = 26% Therefore, the signal is detectable 
26% of the time. 
The accuracy in the determination of the percentage of 


time the signal is within any specified frequency range can 
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be increased by making the voltage amplitude ranges smaller 
and shortening the duration of the time intervals. The 

next diagram shows sample values for the same signal, but 
with amplitude ranges of |]/2-volt and the time interval dura- 
tion halved. It is constructed in the same amnner as the 
previous plot. A check to verify its accuracy will help you 


evaluate your understanding of the procedure. 
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As amplitude range and the time interval length are made 


smaller and smaller, the plot approaches a continuous curve. 


continuous curve 
approximation 
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Mmost practical applications a formula for the continuous 
curve can be determined, and the corresponding function inte- 
Meeeea OVer the amplitude range of interest to find the per- 
centage of time the signal was within that range. 

The probability density function of a signal represents 
the “probability” that a signal will be within an amplitude 
range. As the plot constructed above represents |00% of the 
time of the signal, so too does the probability density func- 
tion represent the total probability of the signal. It is 
Scaled from 0 to |, i1.e., the area under the probability den- 
Pmeevetunction curve is equal to ]. If the detector in the 
example just mentioned detects, with certainty, any voltage 
amplitude of 3 volts or more, and the signal spends 26% of 
the time at 3 volts or more, then the probability that the 
Signal will be detected in the 3-volt or more range is 0.26. 
The probability density is expressed by a function which may 
be integrated to find the probability that the signal is in 
any prescribed amplitude range. For probability density func- 
tions used in many applied problems, the integration process 
can be avoided, sincd the results are readily available in 
standard tables. 

The probability density function (pdf) most commonly used 
fer electrical signals is the "gaussian" or “normal” pdf. It 


is the familiar bell-shaped curve described by the equation 
- (x-m) * 
i} De 
BR Gl 2h 





fiaecie equation for p(x), 


. | ee 
m = mean, or average, value, ™ 2 
15 


x value of the signal, x(t), at time 


a 


fo) standard deviation (o is the Greek lower case sigma). 
The standard deviation, 0, 1S a measure of the spread of the 
Signal abour its average value, m. In the next figure we 


we show two gaussian curves with the same mean, but different 


values of o. The area under each curve is equal to |.0. 
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Por the gaussian pdf, about 68% of the area under the curve 
lies between m - o and m + o, which means that the probability 
that the random signal, x(t), assumes a value between m -lo 


and m+ jlo is 0.68. 
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m+ 2c, P(x) = 0.95 Le 
(Shaded area = 95% ‘CY 
cA y, ra 
of total area) es nes 
we vn ~~ oo 


Bly 





Pe 30, P(x) = 0.997 
[Saaded area = 99.73% 


of total area) 





In the illustrations on the previous page, only 0.997 of 


the total probability has been accounted for. This is because 


there are no limits on a gaussian pdf; the range is unbounded. 
A signal may conceivably exist at any finite amplitude above 
or below the average value. The probability that the signal 
will occur outside the m zs 20 Fange tseonly 0.003 (|). 0520997 
= 0.003). Limiting circuits are normally used early in the 
orocesSing stage to prevent the rare occurrence of a large 
amplitude signal from damaging equipment. The important part 
Siete signal, 0.997 of it, is still retained. 

The gaussian pdf is centered around an average value. 
For an electrical signal this is the average value of the 


amplitude of the signal. Look at the signal and its pdf 


side by side with the pdf oriented to the signal amplitude. 
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The value chosen for m is the statistical average of the 
29 sample values above. Each of the sample values was 0, l, 
Mpes, Or 4. The average, mM, is equal to the sum of the indi- 


vidual sample values divided by the sample size, n. 


Meee ee tte ti tl 24242424+2424242434...4+4... 
29 


or, since many values are the same, 
Inet) + 82) + 7(3) + 4(4) 


—| | oo °° °°.°.°»« « = 
x il nN nN 
Algebraically expressed, this is m= =a al Es where the oa 
means the sum of all the samples, ae PRONG] = eco jo — sna. Or 
xy a Xo iP X2 ata tiecten ct: x: If the number of sample values becomes 


very large, the sum of the products of the individual sample 
values and their probability of occurrence becomes a better 
eeomce for m. Increasing the number of sample values gives a 
Meee accurate description of the signal. If mis calculated 
for the plot with the amplitude range and time duration halved, 
mea=74 and m= 2.095. The algebraic expression, if n is very 


large, 1S approximately 
fete ene. eater) 
Bade ees sec 


Tt has been shown how statistical characteristics may de- 
emibe a Signal. The next problem is how to obtain signal 
Smeracteristics that are statistically valid to identify a 
eweonal . 

The validity of statistical methods in determining the 
outcome of a process depends on repition of the process over 
and over saan. With a random signal, x(t), this repetition 


may be accomplished by using an "ensemble" of random signal 
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generators with each signal generated having the same broad 
meaeesotical characteristics. If the individual outputs from 
this “ensemble” are measured simultaneously at some time, t, 
and averaged, tne result would be an "ensemble average", in- 
dicated by x(t). This "ensemble average" is identical to, 
and has all the properties of, a statistical average. How- 
ever, most signals of interest in an ASW environment originate 
from a single source. Therefore, only the time average of 
the signal, Kx(t) > , May be obtained. 

mihe time average, &(t), andetnevcnsen> lemaverage, x(t), 
are not always the same. For example, suppose the statisti- 
cal characteristics of the signal from the random signal gen- 
erators in the ensemble vary with time. Such a variation 
would not be reflected in Measurements Made at a fixeca time, 


ie The ensemble average at time ty would be different from 


1° 


the ensemble average at time t A random signal with these 


5° 
Seeeacteristics is called "non-stationary". If a detector is 
used to detect the random signal in our example above, what 
Mmemptne effect on the detector if the random signal is “non- 
Stationary"? The average, m, was calculated as 2.]4. Suppose 
that m varies with time, and that over another time interval 
a is calculated to be ].]4. If the signal over the second 
time has the same characteristics as over the first, the 
effect on detection would be the same as raising the thresh- 


mia to 4 volts over the first interval, as indicated in the 


next diagram. 
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The probability that x = 3 volts or more would become 2/29 = 
0.069 vice 0.26 for m = 2.]4! If the average is unchanged 
and the standard deviation increased by a factor of three, 
the detector would be observing only 0.68 of the signal as 
opposed to 0.997. Therefore, if the random signal is not 
Beeationary”" (i.e., its statistical characteristics do not 
vary with time) or at least "stationary" for the processing 


time of the detector, the detector will be very inefficient. 


Even though the signal may be stationary, the time and ensem- 


ble averages may still be unequal. 
Fortunately, many random signals that mjst be processed 
in ASW have time and ensemble averages which are identical, 


perhaps not for all time, but at least for the processing 


time required. Signals with their time and ensemble averages 
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equal, <x (t)> = x(t), are said to come from an "ergodic pro- 
mes ,Or to have the property of “ergodicity". An "ergodic 
process" is also stationary. Under the assumption that a 
Signal from an ergodic process means the time and ensemble 
averages are equal, the signal is stationary and the time 
average will have all the properties of a statistical average. 
” Having extablished that a time average will be statisti- 
cally valid if it is from an ergodic process, how do we ob- 
tain the other statistical characteristics which describe a 
Signal? All acoustic detectors convert the signal (and, un- 
fortunately, the noise, which will be addressed later) into 


some time-varying voltage or current, x(t). 
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The time average of the signal, &&(t)> , assuming ergodi- 
city, is equal to the average value, m. This average value 
is both the value about which the signal's gaussian pdf is 
centered and, as the time average, represents the "dc" compo- 
Meme of the signal (direct current, also called the steady- 
state or the non-time-varying component). If m= 0, there 
768 no "dc" component. If the signal, x(t), is impressed a- 
cross a one-ohm resistor, the instantaneous power dissipated 
mex (t). The time average, GPE) Y, becomes the "second 


moment", which represents the total average power, x°, of the 
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electrical signal. The “second moment" differs from the 


Z 


"average value squared", m°, which is the dc power, or the 


power in the non-time-varying component: m* is obtained by 


squaring after averaging, whereas x* is obtained by squaring 
before averaging. 
As an illustraion of these values, consider a random 


Begmal, x(t). 
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The following drawings illustrate computation of m* and 


* 


MO) 


x (¢) 





There sre two other measures which identify a specific 
gaussian pdf besides the average value, m, about which it is 
centered. A random signal is composed of a dc component and 


it te 


faeeac §6<COMponent (alternating current, or the time varying 


Semoonent of the signal). The total power, x*, must be the 
sum of the dc power, m*, and the ac power, o*. Though repre- 
Memeimg the ac power of the signal, o* is the "variance" of 


the gaussian pdf. The variance is equal to the second moment 


[Mus the average squared, o* = x* - m*. For an electrical 


Signal the variance is the total average power minus the dc 
power. As the standard deviation, 0, represents the square 
root of the ac power, it 1S more commonly known as the root- 
mean square (rms) value. We summarize these definitions 
briefly. 

Given an ergodic random signal, x(t): 


}. the average (mean) value, m, 1s its de component; 
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2. the average-squared, m*, is its dc power; 


3. the second moment, x*, is its total average power; 


feeethe variance, o* = x* - m*, is its ac power; 


Bee tme standard deviation, 0, iS its rms value. 

mi Fourier Transform Properties", the concept of correla- 
tion of a signal was discussed. It is considered again here 
because of its unique properties, and because of its relation- 
ship to the "power spectral density” of a random signal. 

Those who Feel they have not grasped the concept of correla- 
tion should review that section before proceeding. 

The autocorrelation function, R(t), was defined previously 
ia@eeehe time average, R(t) = Cee: + ey For an ergodic 


random signal the autocorrelation function becomes: 


R(t) = C(t) » Sate ot t)) = x(@) X(t <P a) 
Mae value of R(t) is a function of the time delay, T. At T=0, 


R(t) is a manimun: 


R(O) = x(t)+ x(t+0) = x* = total average power. 
As mirereases, R(t) decreases, until at T= infinity, R(T) 


is a minimum: 
R(2) = m-> = de power; 
R(O) - R(e) = 0% = x* - m° = ac Nowe: 

The autocorrelation function is also a measure of the 
"time coherence" of the random signal. If t is small, x(t) 
and x(t+t) will be close together in time. Therefore, the 
presence of one will have some effect on the presence of the 


other, a condition known as "statistical dependence". If the 
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presence of x(t) has no effect on x(t+tT), such as happens if 
tT is large and x(t) is random and non-periodic, then x(t) and 
x(t+1T) are said to be “statistically independent". If fora 


Signal with a gaussian pdf the absolute value of the auto- 


correlation function minus the dc power equals zero, | R(x) ss 
m? = 0, then x(t) and x(t+t) are statistically independent 
and "uncorrelated". Statistical independence simplifies the 


calculation of the autocorrelation function. Whether or not 
two signals are "correlated" becomes important in the calcula- 
tion of their average powers. 

Assume that a signal has been formed by the addition of 
meyenseparate signals: z(t) = x(t) + y(t). The correlation 


mre cion Of z(t) will be of the form 


(xty) (x+y) = x? + xy + yx + y*, and is equal to: 
R(t) = €x(t)+ x(t+t)> + &y (t): y (ttt) + Kelty (ett + Y(t): x(t+t) 
R, (t) ~ ae, tf og Ee a R(T) a peut! 


where the last two terms are the "cross correlation terms" of 
faeeyeewith y(t) and y(t) with x(t). As noted previously, the 
above correlation functions represent the average power in 
the signals, and can be written as: 
oe + Py ar Pay af Bx” 
if the two signals, x(t) and y(t), are correlated over the 
time delay, T, then the total signal power, 5, Toe e Ss Unsor 


the signal power in signal x(t), S.., the signal power in sig- 


XK 


mal y(t), Sor and the cross correlation of their combined 


Signal powers, = ay and Sop! 
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Zz" xX yy XY ys 
Pietne Signal 1s the sum of mutually uncorrelated signals, 
aS in the case of noise uncorrelated over the time delay, 1, 
then the total noise power, N, is the sum of only the noise 
power Ny Miele tne Cross correlatuvons are equal to zero: 


x BS 
N =N +N; (N = (a. jn = 01. 


Z. xX VY XY yx 
For a signal which is the sum of mutually uncorrelated 
meats, the correlation function is the sum of only the auto- 
correlations (the cross correlations equal zero), and the av- 


erage power is the sum of only the average powers of the 


individual signals (the cross correlated powers equal zero). 


EeeeeeOWER SPECTRAL DENSITY 


Mos fai only the time dao 
been addressed. The frequency domain is represented by the 
power spectral density, which was introduced in "Filters and 
Bandwidth Requirements". The device used to transfer between 
the time and frequency domains is the Fourier Transform. The 
Weiner-Kinchine Theorem states that the power spectral density, 


mee, and the autocorrelation function, R(t), are Fourier 


mieanstorms of each other: 


f- 00 


G(f) = FiR(t)] = f x(t) x(t+re 


a eelae 


And, from the Duality Theoren, 
ieee (1) <>Git), then G(t)+> R(E£). 
As can be seen above, the random signal may be described 


in either the time domain,’ R(t), or the frequency domain, G(f), 





fmt tree interchange between the two. 

Since the power spectral density is the Fourier Transform 
Seene autocorrelation function, the Fourier Transform Theo- 
rems may be applied to the power spectral density. The sys- 
tems designer uses the properties expressed by these Peotone 
to design an efficient detector/signal processing system. The 
ASW specialist who has knowledge of these properties can bet- 
ter understand how his system is designed to operate and 
therefore, utilize it more effectively. The duality theorem 
was already utilized. The other major theorems applicable to 
power spectral density will be explained along with examples 
of their usage. 

ie) 6©Frequency Translation Theorem 

ieee sicemal, x(t), 1s bandlimited in W < ee and a 
mewesiqnal, y(t), is formed by multiplying x(t) by 

cos (wu tt>), Watt) = x (e) cos (Ww tt¢), then the correlation 

function of y(t) is equal to one-half the correlation 

miaction of x(t) COSWT, thateis Ru = (]./2)R, (1) cosw ot. 

In addition, the power spectral density, G(f), is 


G(£) = (1/4)G, (££) + (1/406, (E+E) 


Multiplication in the time domain by cosu) t becomes transla- 

tion in the frequency domain, as the power spectrum is shifted 
up and down in frequency by +f. THemecone ii TOnerenat tiie jatace 
width, W, be less than f is necessary to ensure that G(f£-£ 0) 
and ee itt 7) do not overlap, otherwise aliasing, as described 


in "Filters and Bandwidth Requirements’, would result. Notice 


oS 





feo that in multiplication by cos (wu ttd), tne phase factor, 
@, 1s lost in translation to the power spectral density. 
Mi@vefore, the original random waveform may not be recon- 
structed from knowledge of the power spectral density. How- 
ever, by returning to the time domain via the inverse Fourier 
meanistOrm, the autocorrelation function, R(t), with its sig- 
miertcent Statistical properties, may be found. 

Frequency translation is very useful in signal processing. 
Suppose that due to space, weight, or power limitations it is 
impractical to locate the detector and the processor together. 
The low frequency (0 to J000 Hz) sound may be detected, trans- 
Yated to radio frequency (RF) (HF = 2 - 30 MHz, VHF = j]00 - 
200 MHz, UHF = 200 - 400 MHz) and broadcast to the processor. 
In the same manner, one processor can handle many detectors 
if each detector's signal is translated to a separate fre- 
quency band within the RF range (Detector |, frequency fy loko) 


if Detector 2, frequency f., Ee. fyi Gibe wal) wathane toe sUHE 


2! 
band). The problem could be one of processing signals of 
many different frequencies. Rather than having many differ- 
ent processors, each processing only one signal, it may be 
Metter to have one processor Operating at a fixed frequency. 
mot the signals could be translated to the processor fre- 
Guency for processing. If the processor 1s very elaborate 


and/or very expensive, this method becomes desirable. 


2. Integration Theorem 


t 
= ' ty il fe 
ety (tt) -f Mit wicitne hen Eye) ae G,, | ie 
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The multiplication factor, (anf)? ? indicates that the low 
frequency components of the power spectral density will be 
enhanced by integration. The ocean acts as an "integrating 
filter", in that it allows the low frequencies to propagate, 
while rapidly attenuating the higher frequencies. This is one 
reason why, in an attempt to get long detection ranges, active 
Semeaes and passive processing systems have increasingly ex- 
ploited the lower frequencies. | 

3. Differentiation Theorem 


This theorem is offered without an example since it 


follows so naturally the integration theorem. 


ire e— d[x{t)}/adt, then G,, (£) = (2u£) °G. (£). 
Geeemultiplication factor, (Dee, sbenvcbleeteee ete ae lnaeya 
frequency components of the power spectral density will be 


enhanced by differentiation. 


C. NOISE 

At the beginning of this section, noise was mentioned as 
a "special case of random signal which is always present to 
some degree whether or not a ‘target’ is present". In com- 
munications, noise is often defined as any electrical inter- 
ference, or unwanted signal. Signal then, in aAdgveren) tO the 
characteristics previously used to describe it, has the qual- 
ity of being wanted, sought after, emanating from a source 
of interest, and containing intelligence. In ASW, signal is 


always associated with "target". Noise, on the other hand, 


at 





2S anytning which interferes with, or tends to. mask, the 
megial. Often “one man's signal is another man's noise" as 
any patrol aviation Tactical Coordinator whose buoy field has 
miisiasoeen penetrated by a "pinging" Destroyer will verify. 

An acoustic detector converts both signal and noise into a 
time-varying electrical voltage or current. An examination 
of some of the characteristice of noise wiil help in under- 
standing how systems may be designed to reduce it. 

Another way to define noise is simply as that guantity 
observed in the absence of signal. Noisé may come from a var- 
1ety of manmade or naturally occurring sources. Systems are 
designed to eliminate as much noise as possible, but some 
noise will inevitably remain. The thermal motion of electrons 
in the conducting media ot the compvonents of a system, for 
instance, is one unavoidable cause of electron noise. Ther- 
mal noise iS interesting in that its power Spectrum is con- 
stant over a wide range of frequencies. It is designated 
Byaete noise” by analogy to white light, as all frequency com- 
ponents are present in equal amplitudes. The amplitude of 
thermal noise has been proven to have a "gaussian" pdf. It is 
therefore referred to as “gaussian white noise" ("gaussian" 
Mmmpretude distribution; “white” frequency distribution). If 
"gaussian white noise" is played over a loudspeaker, it sounds 
dull and monotonous, somewhat like a waterfall. The subtle-' 
ties of its random variations are hidden from the human ob- 
server. 

The power spectrum of gaussian white noise, with zero 
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: lee ; 
mean, is G(f) = 5 » where n is power density in watts pew 


Hertz. By Fourier transformation, the autocorrelation func- 
tion is < 
R(T) =f ail’ Shae 
~OO 


the "delta function" introduced previously. 


a | 
me | — 


ca ty, 
ess | 


oe , Where 6 represents 


a 


fers apparent from the R(t) plot that for any time delay, T, 
Sener than zero, the autocorrelation function is zero; there- 
fore, any two different samples of a gaussian white noise sig- 
nal are uncorrelated and statistically independent. 

The signals of interest in ASW are frequency limited, 
1.e., there is some finite frequency band within which the 
emomadl exists. Filtering in order to improve detection of 
that signal was discussed earlier. If aqaussian white noise 
is filtered, its frequency components are naturally those of 
the filter bandwidth, but the amplitude distribution remasins 
gaussian. The output power spectrum of white gaussian noise, 


filtered by an ideal filter, is a rectangular function, 


a n 
ot) i (—) (SR) where Il indicates a rec- 


tangular function and B represents the filter bandpass. The 
autocorrelation function is the inverse Fourier Transform of 
Mme power spectral density. As indicated in "Fourier Trans- 


manm Properties", the Fourier Transform of a rectangular 


a7 





-@ | +8 





white (Gt usscow idecl bandpass 
Nocse Ecller 
@mese 1S a "SINC" function, Ba SU eh ouNe2 bt 





The above figure shows that filtering has produced the follow- 
amo results. 


1. The power spectrum, though no longer white, is con- 


Seant Over the finite frequency range of the filter. 


Peeeerne Gutput power is finite, y* = / (2B) = Noe Le 


varies linearly with bandwidth, B. 

fee the output signal is correlated over time intervals 

@f about x. 

How is this knowledge used in signal processing? 

Item (]). The noise in which a signal of interest must 
be detected is often not gaussian white noise. However, 
through judicious choice of a bandwidth, its power spectral 
density may be constant over that range. Otherwise, a pre- 
whitening filter may be used to make the noise power spectral 


density constant over the bandwidth. Therefore, after filter- 


ing it may be treated as gaussian white noise. 
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Item (2). The noise power after filtering is a= Be 
Miers a direct, linear function of bandwidth, B. By narrow- 
ing the bandwidth, the amount of noise power present is re- 
duced. Therefore, the signal power is enhanced in relation to 
Ee noise power. This is crucial in a power, or energy | 
detector. 

Item (3). The bandwidth of the filter determines the min- 


imum time delay, T, in a correlation detector. It must be 


plas ne 
2B! 2B 


fmm ne correlated. If the noise is uncorrelated, which it 


greater than (= ), or the noise, as well as the signal, 


se, only its average power will add in a corre- 


lation detector, thus decreasing in relation to the signal 


Mari be for T> 


power. 
ime bandwidth of the filter also determines the finite 
integration time to be used in determining the time average. 


The time average is indicated by €}(brackets) and is defined 


as: 1 Ah 
lim 2 af x(t) x(t+t) dt. 
Tx 0 
Replacing T by d eyatcl — lah Sa we obtain 
2B ] 
2B 
BJ x(t) x(tt55) dt. 


0 


A review of the development of bandwidth versus integration 
time in the "Filters and Bandwidth Requirements" section may 
be helpful in the understanding of these relationships. 

This section has been an attempt to describe the concepts 
of random signals, power spectral density, and noise. The 


intent has been to provide a basis for understanding their 
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individual characteristics. A grasp'of these characteris- 
tics 1s necessary to understand the methods of Signal process- 
ing. The most common methods in current use, DELTIC, energy 
detectors, correlation detectors and beam formers, will be 


presented in following sections. 


Oe 
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Miethe section on filters and bandwidth, the relation- 
ship between integration time and frequency resolution 
was discussed. It was shown that one must integrate for 
longer periods of time as the frequency resolution 
requirements become stricter. In many cases, however, 
the integration time required to give the necessary 
meoolmtion 1S much too long to be of tactical use. In 
O@rder to circumvent this, several processing schemes 
have been developed which in effect give more resolution 
for a given integration time. One of the more common 
in use with military processors is known as the DELTIC 
(DELay TIme Compressor). Let us investigate the manner 
mewiateh this processor accomplishes this feat. 

it has been shown that processing the signal as it 
imemneecived by the detector will not give any greater 
mesolution. The DELTIC takes the signal as received, 
omeerecords it on a magnetic tape, disc, or drum. It 
fees this in a special sequence, however. To see this, 
let us take a simple signal and observe how it is processed. 


AMPLITUDE 





For simplicity in showing the operation, an 8-sample 
point "time window" will be used. In actual practice, 
1000 or more sample points are commonly used. If the 
time between sample points is %, the time window observed 
feeeeeoe St = T. Therefore if the DELTIC is not used, 
the best resolution we can get is 1/8%T, and it would 
mauere St seconds to get this resolution. The physical 
layout of a DELTIC is shown in the next diagram. 


INPUT READ: GUT 


“GATE” , _ L HEAD 





Secs? 


11 S\N A. MAGMETIC DRUM, 
a DISC, ER TAPE 

iite drum is rotated at a rate such that it makes one 
MevoOlution plus one sample space during the time period 7. 
eens case, with 8 sample points it makes 1.125 rev- 
Meterens per ~t seconds. The read-out head continuously 
Meads OUL the sample values as the drum rotates under it. 
fie input gate injects a new sample point once each 
meres revolutions in such a way that the sample point 
which has been on the drum the longest time is replaced. 
Mearting with a blank drum, let us examine the process 


es it occurs. 


1 Oly 








READ-OVT BURiIMG FIRST 
REVOLUTION 


_— > | - Beso 






When the drum has made 1.125 revolutions (in SeCG.4 


tne input gate inserts sample point #2. 


SAMPLE 4a 








peeeer 1.125 more revolutions, the 3rd sample point 


ime inserted. 








eile remaining spots are filled similarly until all 


eewestilled. The output during one revolution of the 


jete[ste[s[elz[e] 


mecimeene next Sample point is inserted after 1.125 revo- 


arum then is: 


mimeons, it displaces the first. 





Once the drum is filled (which takes 8% sec.), a new 
foe replaces the oldest one every tT seconds. If we 
HOok at the output from the read-out head, we see that 
there are 9 values read out each t seconds (8 values/ 


revolution, 1.125 revolutions/t seconds). 


a Set. — 


Mmiemwoutput of the DELTIC each t seconds is thus the 
Same as the inputs over the last 8% seconds. The DELTIC 
operation described above can be illustrated by the 
following analogy. I1f we record 80 seconds of conversation 
on tape, and then speed up the playback so that the 
conversation takes only 10 seconds, we have compressed 


Mmenconversation just as a DELTIC would. The same 
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meormation 1S in the 10 second playback as was in the 
80 second conversation, but it only takes 1/8th the time 
to hear ieee also notice that this compressed conver- 
Sation has been shifted in frequency by the compression, 
and normal voices sound like Donald Duck. This illustrates 
@emeractical effects of the Fourier scale change theorem. 
Beetising the DELTIC we then are able to observe the 
preceding 8t seconds worth of signal during each interval 
of length t seconds. The frequency domain plots of the 
Original signal and the output of the DELTIC show that 
mivewOriginal spectrum has been shifted in frequency and 


Semetched out. 


| CRIGINAL SPECTRUM 






Minteium RESOLVABLE 
“7 FRCOvaNcy BreFCRENCE 


DELTIC OUTPUT SPECTRUM 








The Fourier scale change theorem shows that we have 
foeeerplited the time domain characteristics of the signal 
By the factor 1/9, and so the frequency domain character- 
meres Of the signal have been divided by this factor, 
Mao ecrtitectively multiplying the frequency plot by 9. 

iiemecompression factor" is the ratio of the original 
aampee time, [, to the compressed sample time, T'. 
imo in this case it is the ratio T to T/9, or 9. Note 
meaeenot Only is the center frequency and bandwidth 
multiplied by this factor, but also the minimum resolvable 
frequency difference. Without compression, we would 
me oatate tor T seconds to obtain a resolvable difference 
of 1/T. Upon compression, this becomes 9/T. At first 
feeeece, this does not seem to help, since the resolution 
feeenot been improved. Why use the DELTIC if it doesn't 
nomep ¢ 

The benefit in using the DELTIC is that we can obtain 
this resolution not by integrating over I seconds, but 
mer) 1/9 seconds once the DELTIC has been filled. Many 
em@ectrum analyzers ‘'scan" the frequency band of interest 
byeusing a single band-pass filter of bandwidth equal 
@omthe minimum resolvable frequency difference of the 
Mageessor, and stepping its center frequency across the 
band one step each integration period. Thus, if the 
filter bandwidth is 1/9 the pameiden of interest, without 
the DELTIC it would require 9 integration periods (T) for 


one scan of the band. With the DELTIC, each integration 
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Meriod 1s 1/9, and thus we can scan the entire band with 
the same resolution in what was previously one integration 
Meme, !. This 1S quite an advantage for a band of 
Several hundred hertz and a required resolution of 0.5 Hz. 

iirc DELTIC need not be designed for insertion of a 
fevesample value every revolution, but only a new value 
moproximately every 10 revolutions. In this way one can 
Bompress the signal even further by having the drum 
mevolve more than once every %?% seconds. 

Instead of a magnetic drum, disc, or tape memory, 
meceeital Shift register may be used. The shift register 
Mmepsecppcd past the read-out gate in such a way as to 
read the entire register each t seconds. The values 
meemwenen recirculated. Each time a complete circuit is 
meeee a new value is input to replace the oldest, and 


mmemprocess continues. 


RBAGEE 
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meme systems which use the DELTIC tn processing are: 
a. AN/AQA-7 
bee AN/ BOR-1 5 
ee iN / SOS- 26 
d. AN/BQR- 20 
Meet litustrations on the following pages show a side hy 
Meeemeconmparison of an input signal and the resultant DELTIC 
Meee NOtiCe the gradual build up of the DELTIC output as 
Megemeincervals of the signal are introduced. Note also that 
Smee intervals continue to be introduced, the DELTIC drops 
Seems Oldest information and represents only the most re- 


mre cime intervals. 
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fPemereesional 15 “Stepped” into the DELTIC, the output 


—- 


eet 


Seeene DELTIC represents the compressed signal. Plots A to 
memaieate the buildup of the DELTIC. Plots G to J show the 
eldest information dropping out and only the most recent in- 
formation being represented. Plot K shows the complete time 


femeery Of the signal and its DELTIC output. 
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In the discussion of random signals it was noted 
miae in Order to apply statistical analysis to signal 
processing, the signals must be "ergodic". One property 
meemereeodic Signals is statistical stationarity, that is, 
miemstatistical properties of the signal are constant 
miecime. In reality, the signals that must be processed 
are not strictly stationary. The sea noise varies with 
mmemcime Of day, changes in weather, etc., but the var- 
moron iS over a long time period. This fact enables 


memeto apply the statistical analysis approach to signal 


« 
Cy in te a ae re 


Mee ssing, aS long as the proccssing time is short cncug! 
mee the Variation in signal properties is negligible. 
Let us consider how to apply statistical methods to 
detection problems. For simplicity, take the case of 
Gaussian noise, where the probability distributions of 
amplitude for noise alone and noise plus signal are 


: . : Z : 
Mmissian with equal variance, 6 , but different means 


(average values). 


~ 
4 





If we now have a processor which compares the amplitude 
Seeche input to a "threshold" level set into it, and 
maidicates “signal present" if the input is above threshold, 
“signal absent" if the input is below threshold, then 
we can apply statistical theory to determine the oper- 


femme characteristics of the processor. 






we q 
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THRESHOLD 
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The area under the S + N curve above the threshold 


+ 





Boeeproportional to the probability that the amplitude of 
the input will be above the threshold when the signal 
mamemescnit. tIhus, this area is proportional to the 


meena! detection probability, P(d). 








The area under the N curve and above the threshold 
is proportional LO tmcmpLOvAba lity that) the amplitude: of 
the input will be above the threshold, even when the 
Signal is absent, because of the background noise. Thus, 
mums area is proportional to the probability of false 


feearm, P(fa). 





*£ (cs) 
ff the threshold is varied, the probabilities change. 
Meee threshold is raised, the probability of detection 


Mereases, but so does the probability of false alarm. 





Pemtne threshold is lowered, P({d) increases, but 


samdoes P(fa). 
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Mieeratio of P(d) to P(fa) does not remain constant, 


however. Plotting P(d) versus P{(fa) as the threshold 


Varies, we obtain the foliowing Tresuit. 





E ee 


iimuseis known as a ‘Receiver Operating Characteristics" 


mame or ROC curve. Using this curve we can determine 


Milere to set the receiver threshold to optimize P(d) 
and P(fa) 
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As the threshoid is reised, P(d) .does not change 
meemein Comparison to the rate at which P(fa) is 
decreasing, until the "knee" of the curve (between point 2 


mace point 3) is reached. 
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Beewe continue to raise the threshoid above this 
fee, P(d) starts to fall off rapidly, while P(fa) 
changes very little. By deciding what maximum P(fa) 
can be tolerated, we can set the threshold to yield the 
maximum P(d) possible for that particular P(fa) by 
jewetring to the ROC curve. Conversely, if a minimum 
P(d) is required, the ROC curve indicates the minimum 
P( fa) that must be tolerated. 

Until now the effect of changing the mean of the 
Seren distribution with respect to the N distribution 
would have on P(d) and P(fa) has not been discussed. 
imemtne mean of the S + N distribution is raised, then 


the waren Eade P( fa ) icmlaneer tor alimtinmeshold settangs. 
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The result is shown on the ROC curve in the following 
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ee 
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diagram. 
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Wemdertine the "detection index" as the quantity 


GX - ie (6, which is a measure of the difference 


Sth 
between the distributions. A large d value indicates a 


large separation between the distributions, and thus a 


Sieve with more of a ''knee''. 
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P(d) 





Para) 


mais diagram shows that if the presence of a signal 
etaeically changes the input distribution mean, then 
Memean reduce the P(fa) Seon camel ewan lews ta. lel 
em cmue ) Cae Ii whe oresence cf a: sicnad 
does not drastically change the mean, large values for 
P{d) are possible only at the expense of a large p(fa) °. 

We have discussed thresholds, and the way the selection 
Mememmeshold values affects the probability of detection 
migeeene probability of false alarm, but have not discussed 
Mewethnese thresholds are set. One point that may avert 
bemreadily apparent is that in detection systems in which 
Giemnuman operator is involved in the decision as to 
Presence or absence of a signal, his perception of the 
ilay must be included in computation of the threshold 
value. Most systems in use eee oneal, PiweOLVe din 
Operator in the decision process. The SQS-2Z3 sonar is 


an example. The operator can vary the threshold to an 
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extent by adjusting the brightness and intensity of the 

display id the gain control, but a large consideration 

in determining whether or not a Ba Weetels present 15) his 
ability to distinguish between target and noise as pre- 

memeed On the PPI scope. 

The operator reading grams on an AQA-7 is another 
eeeuple. He can control the gain, or the grey value of 
tne background on the paper, but his perception of the 
distinction between target and noise determines the 
threshold. For this reason it is Qi feu 16 sie we muMlys 
the ROC curves for systems with human operators, because 
iteas difficult to determine the actual threshold setting. 
The threshold will vary from day to day, or even hour 
to hour, with the same operator and the same hardwd:ie 
Settings because of the operator's state of mind, physical 
Gendition, fatigue, etc. This is what makes quantification 
Seepeie detection characteristics for a system difficult. 
Samioletely automatic systems are much easier to evaluate, 
mm—ome can determine exactly what the settings are. 

Im the preceding discussion of energy detection and 
thresholding we have seen how a system's detection capability 
can be determined from knowledge of the characteristic 
@istributions of noise and signal plus noise. Implicit 
in this treatment is the assumption that the distributions 
of noise and signal plus noise are exactly known. If these 
oistributions are exactly known we can specify the per- 


formance of the system for any given threshold setting 
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mmeterms of P(d) and P(fa)'. If we do not know these 
Meter? butions exactly, our analysis will not exactly 
mreenlbe the system performance. Thus, the accuracy of 

our prediction of system performance depends on the 
accuracy with which we can describe these distributions. 
paec@esipning real world systems this is one of the problems 
facing the designer,since the distributions of noise and 
Signal plus noise are not constant in time (over any 
memetiy period) or with geographical position. For this 
reason many systems do not always perform as well as 
memespecifications indicate. Only if the ambient conditions 
are the same as those for which the system was designed 


feet perform as predicted. 
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rie CORRELATION DETECTION 


In the section on random signals, power special 
density, and noise, it was shown that correlation can 
er oemiyety improve the signal-to-noise ratio for cor- 
meeerea Sipnals due to the cross product terms in the 
meeemelation Of the correlation function. How can this 
marormation be used to detect signals in noise? Let 
Memeceudy an active echo ranging system as an example. 

In this type of system a signal is transmitted, and 
then the returning signals are processed. The range 
Of the reflecting object can be determined from the echo 
return mimes sino a noise-free environment, there are 
few problems in this system, but in most real situations 
micwambient noise may ‘bury" some of the returns. By 
Using correlation, we can improve the signal-to-noise 
jmeeo ana recover these echoes. 

This system retains a replica of the transmitted 
Signal and then inputs the received signal into the 
correlator. The output of the correlator is large when 
Pamecho 1s received, and small when only noise is present. 

Mine does doppler do to the correlation of returning 
echoes? Let us look at a burst of CW signal transmitted, 
mmgethe effect of doppler shifting on the correlation 
function. If the echo is not doppler shifted, the output 


of the correlator will be the autocorrelation of the 
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transmitted signal at a time , which is related to 


mitemrange of the reflected object. 





Be ellen GOitp tees a Ot 
2EING output 


fameteeneturning echo is shifted in frequency by 
moppler, the output of the correlator becomes the cross- 
Semeclation of two closely related signals, but its value 
1s never as high as the autocorrelation. The cross- 
Semrelation amplitude and the amount of delay over which 
mammprcmcorrelated decrease. (In fact, the reduction results 
feme uecrease in the delay over which the Sipnal is 


eempelated with its doppler shifted echo.) 


i | 


iEanscmi t ted doppler shifted EO wits lle eee 
Signal returning echo output 





Pepletsetmtne Correlator output for a number of different 
aeopber shiits on the same diagram, that is, av, Af, 


Seerelator amplitude plot, shows the effects of doppler. 
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meme Given echo, this®plot gives a solid surface 


relating Af, %, and amplitude. 
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EQuAL AMPLITUBE 
CONTOUR LINES 


ieewe look down at the ‘T,Af-plane, and plot the 
contour line where the amplitude of the function is 
Bumper cent of the maximum, we obtain the following 


result. 
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This diagram is commonly known as an "ambiguity 
diagram", and the reason for this will become apparent 
Peemeexamination of FM pulses. First, however, let us 
consider a long CW pulse to see how doppler affects the 
merrelator output. 

mmeesection On correlation and convolution shows 
Mememcie width of the correlation function is related 


Momence duration of the time-domain signal. 


AUTOELCORRE LATION 


TIME DOMAIN | FUNG TION 
SIGHAL 
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Maus, £Or a iong CW puise the correiator can oe 
Beeeeected £0 respond over a long range of values of 
eieanout doppler and its effect on a long CW pulse 
being correlated? Observation of two signals at slightly 
different frequencies shows that the longer the pulses 
memme correlated, the less the correlation, since tne 
emails Will be almost correlated over only a few cycles, 


mimeeet turther and further out of phase as time progresses. 





As the difference in frequency between the signals 


increases, the number of cycles over which the two are 
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meilost Correlated decreases rapidly. Thus for long 
mepuises, the correlator response decreases rapidly 
meemedoppler shifts. The ambiguity diagram for a long 


miepulse is shown in the following figure. 





ime length L is proportional to the length of the 
Giepulse. Tp is the delay corresponding to the actual 
fave of the target. 

eee Short CW pulse, the length of the correlation 
medio) per cent contour will be shorter, but the range 
Semeoppler shifts over which the signals remain essentially 


@ermclated is greater. 





| Purse 
LENGTH 


ims the ambiguity diagram for the short CW pulse 
has the following appearance. 
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These two diagrams show why long pulse CW is good 
femeesecarching and doppler determination, and short pulse 
meets) fO0d for range resolution. A sonar pulse of 120 
milliseconds will ensonify a region of water 180 meters 
iomee in other words, the spatial length of the pulse 
as it traveis through the water is i806 teters (1500 meters/ 
Second velocity x 0.120 seconds = 180 meters, the distance 
the forward edge of the pulse has traveled during the time 
the rest of the pulse was being generated). The cor- 
fmemrtion peak 50 per cent contour, as shown on the 
ambiguity diagram, will be proportional to the pulse 
memeth T, and for CW this turns out to be a one-to-one 
mitro. Thus, the correlation will be over a period of 
iwemilliseconds, which means that the target range can 
be determined only to the nearest 180 meters. 

Rmmememeise ot the siiort pulse of, say S milliseconds, 
the peilccespatial length is 7.5 meters, the Po ciatien 
Medks over 5 milliseconds, and the range resolution is 


fe> meters. 


127 





The effect of doppler on the ambiguity diagram shows 
that the same type of development can be made to relate 
frequency resolution of the two different pulse lengths. 
Because the long pulse remains correlated over very 
limited doppler shifts, the frequency of the returning 
echo can be determined rather accurately by using a bank 
of correlators comparing the echo to replicas with varying 


eounecs Of doppler shifts. 
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On the other hand, use of the same bank of correlators 
feeemea short pulse produces a response from each of these 
Semrelators, because the short pulse remains correlated 


emer a large range of doppler shifts. 
af | 
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If we use a pulse in which the frequency is changed 
during transmission, which is commonly referred to as 


BMeslice, several differences appear. 


Ag aunetegqrntny 
FM "Up-Slide" Waitt 
| | | 


ai 


i avis ull 


Because of the frequency change during the pulse, 
Miemreturning echo correlates with the transmitted replica 
Mer a Short range of time, €, even for long pulses. 

Thus, the ambiguity diagram is very narrow in the % 
MMseti0n. On the other hand, if a ionug puise is used, 
the Af range of correlation is narrow. One difference 
between FM and CW pulses becomes apparent here, however. 
meine Case of FM pulses, a doppler shift will tend to 
amet the entire pulse. in frequency, but the slide will 
remain linear (provided it was originally linear), and 
mits at Will still tend to correlate but at a different 
welue of tT. This behavior and its effect on the ambiguity 


diagram are shown in the next drawing. 
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Note that the FM slide pulse has some of the char- 


memeristics of both long and short pulse CW. At any 
Sven doppler shift, the correlation peak itself is 
iemeow, as 1S the peak for short CW pulseS# Thus, range 


Peaolution 1s good. 
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==WIDTK OF 
GGORRELATION 
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GIVEN DOPPLER 
SWIFT af, 
Note also that at a given vw the range of doppler 
shift over which the signal is correlated is narrow, 
on the order of the range over which the long CW pulse 


[memmeorreiaied. 
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OVER PoePPrer 


ihtesteason this diagram is called an “‘ambiguity 
diagram'' can now be seen. The correlator output, that is 


the v-amplitude plot, shows the following. 


eh 
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Was the peak at To caused by a target exhibiting 
iemaeppler at range equivalent to To» or was it caused 
by a target with 5 knots of doppler at range T, 21 
Miemoutput looks the same in either case. There is a 


Saeomee ambiguity'’ present. 
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One way of resolving this is to transmit two FM slide 
meses, one with an “up-slide", (increasing frequency 
with time), and one with a "“down-slide"(decreasing frequency 


with time).By then using two correlators, one with an 
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mest ide” replica and one with a ''down-slide" replica, 
the ambiguity can be resolved. The next ambiguity diagram 
shows the output from both correlators for a given range 


target. 
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limem compared on a T-ainplitude plot, the correlator 
Outputs for a target exhibiting the doppler shift shown 
on the ambiguity diagram by the dotted line appear as 


in the following diagram. 
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Ey emoting which peak appears first, we can determine 
metne doppler is up or down. If the “up-slide" peak 
memears first, then the returning echo has been shifted 
down in frequency, and thus correlates at a delay, , 
Sieerer than it would if not doppler shifted. This 


midicates down doppler. 


| 
| Ae ue 1M ee "ts 
HTT e st He 
vi i | REPLICA 
| ai 
| ' q tj Hail hiiaiiel ihietall i rity 
h UEAALAPAS Ny 
bila i ee 
We ECHO 





Hl 
Audit Hy 
| 





merthe "down-slide'' peak appears first, the opposite 


[emerue, which indicates up doppler. 
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The true range of the target is midway between the two 
correlation peaks, regardless of the amount of doppler 


mere present. If the target has no doppler, the two 


feanks COINCide. 
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The separation hecueen the up-slide and down-slide 
Somrelacion peaks is proportional to the amount of doppler 
Mmmm the returning echo. If a suitably calibrated 
mroplay of the two outputs is available, the amount of 
Geppler may be read directly from the presentation. 

Thus, by employing FM slide wiles and suitable 
feeocessing, one can obtain the range resolution of short 
CW pulses and the doppler discrimination of long CW pulses 


from a pair of long FM pulses. 
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PweGticr type Of pulse which provides the “best of 
both worlds" with respect to doppler and range is known 
as the "pseudo-random noise pulse", which consists of a 
“coded'' pulse with amplitude or frequency raed a years! 
pseudo-random fashion. Take, for example, the case where 
eemamplitude (or frequency) of the signal is not a 


mole Sinusoid, but a binary coded function. 


AMPLITUDE —-> 
CoR FREQVENCY) 


mis pulse is correlated over a very narrow range 
Meme eiving good range resolution, and it is also cor- 
mobated Over a very narrow range of doppler shift because 
of its pseudo-random character (Truly random signals 
@errclate only at a single value of T.) The ambiguity 
diagram for a pseudo-random pulse is shown in tne next 


figure. 
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Te 


mote that because of the correlation over narrow 
ranges of Af and T, no range ambiguity is present. The 
problem associated with this type of pulse, however, 


= 


is that the presence of noise in the returning ecno 
Causes marked degradation in the correlation of the 
return, and thus degrades the performance of real 
Systems. For this reason the pseudo-random ee USS 
is useful in applications where the signal-to-noise 
feo) 1S high, such as close rafige tracking, bathymetry, 
mee., DUt its usefulness is limited in long range search 
aioe tracking, where the signal-to-noise ratio is much 
Mower . 

PoerOmencmpoint we have covered the use of correlation 
mimaetivye systems, and how the use of FM slide pulses 
combined with the proper correlation processor enables 
iemto determine the amount of doppler in an echo directly. 


Let us now examine some applications of correlation in the 


festive detection field. 
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Saye difference between active and passive correlation 
Memweection Systems becomes apparent very quickly. That 
Mopein the active system we compare the returning echo 
with a replica of the transmitted pulse. In other words, 
miemotonal we are seeking is a known quantity (a replica 
of the transmitted pulse). In a passive system, however, 
iomoulse 1S transmitted, therefore we have no replica 
memecomparison with the incoming signal. Therefore we 
cannot utilize correlation in the same manner; however 
mammerse still] of use, particularly in passive direction- 
mmdane systems. The utility of passive direction-finding 
Systems can be seen by examination of a typical system 
and its operation. Consider a system composed of two 


hydrophones, separated in space by some distance, da. 


od a 


L (A 

If the output of hydrophone 1 is used as a replica 
with which to compare the output of hydrophone 2, the 
@mess-correlation of the two outputs can be obtainéd. 
What does this do for us? The important property of sea 
emcee tor this system is that it tends to be isotropic, 
m@mreis to be the same regardless of the direction from 
which it impinges on the hydrophone. In addition, sea 
noise tends to be uncorrelated over distance. On the 
marer hand, most signals of interest to us are from a 
Single source, usually a submarine, and thus are coherent 


Memomeextent, and directional in nature. 
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pamee SCa nOiSe€ is isotropic, the correlator Output 
will be small and relatively constant when no signal 
is present. Now consider what happens in the presence 
Sieaesipnal of interest. The acoustic TCR 
outward from the source will impinge on the hydrophones 


from some given direction. 


_ 
.-—_.— ++6 
1 2 


Poltowing a given wavefront as it travels across the 
hydrophones, we see that it will encounter one of the 
Myemwophnones before the other unless the direction of 
Mmpenecient is perpendicular to the line joining the two 


hydrophones, or the baseline. 
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iiescaisterence between the times that the wavefront 
mieukes the two hydrophones depends on the direction 
momewhich it arrives with respect to the baseline. If 
miemcanection of arrival is perpendicular to this base- 
line, the wavefront strikes both hydrophones at the same 
mime. and thus there is no time difference. If the 


merection Of arrival is parallel to the baseline, the 
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time difference is a maximum. In fact, the time difference 
Meerclated to the arrivai direction by the following 
equation: 

ce ae 


VecOs 
where v is the speed of sound in the medium, and 6 the 


? 


meevetront arrival angle. 


y 
pr @ 


ite cross-correlation of the two hydrophone outputs 
for a single wave passing the arrav. shows that the delay 
mines t at Which a correlation peak appears is a function 
Seeune Wavefront arrival time difference at the hydro- 
phones. If the wavefront arrival angle is 90°, then 
m@eearrival time difference is zero, and thus the cor- 
memteron peak will appear at T= 0. For any other arrival 
Mimelen there 1S a positive arrival time difference, At, 
paeethe correlaticn peak appears at T= At. his allows 
mS to relate the correlation delay time, %, to the wave 
moma’ direction with respect to the array baseline. 
There is a problem of directional ambiguity, as arrivals 
from opposite sides of the array at equal angles with 


Mespect to the baseline produce the same delay. 
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This ambiguity can be resolved by adding another 
meee tO CLOSS-f1x, or by adding a third hydrophone not 


tiene with the first two, to resolve the ambiguity. 
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It should be noted that although this method works 
well with broadband radiated signals, it has some draw- 
backs when applied to signals of a single frequency line 
or those with a very broad autocorrelation function. 
miesproblem is as follows. If a signal consists of a 
pajabe frequency line, its autocorrelation function is 


@ sinusoidal function. The question then arises as to 


which peak corresponds to the signal arrival direction. 
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For this measOle a COrrelation detector 1s not well suited 
Berepassive detection of narrow-band signals, but works 
meimeeror broadband signal detection. 

Wemiave now seen how correlation can be of help in 
active systems and in passive broadband signal detection 
systems. What can be done about single-frequency line 
etection? In order to auswer this question, we examine 
a system which is similar to correlation in many aspects, 


Mirets NOt actually correlation. 
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IX. BEAM FORMING 


ream torming, or electrical steering of arrays, is 
meeomplished by a process which is similar to correlation, 
but with one major difference--it uses summation, vice 
teciplication, of the outputs. We examine the process 
to see how it operates. 

Let us start with a two-hydrophone array. As noted 
mime preceding section, the acoustic wave arrival direction 
determines the wavefront arrival time difference at the 
hydrophones. With this in mind, suppose that it is 
Meemmea to steer’ the array to receive signals from a 
ME ection rolative tc the arrav baseline selectively, 


mamteis, to form a receiving beam" at a given orientation. 


DIRECTION OF 


ys DESIRED BEAM 


© © 
1 2. 
If we delay the output from hydrophone 2 for a time 


Sige and add t, to the time necessary for the acoustic 


“4 d 
wave from the desired direction to travel from hydrophone 2 
to hydrophone 1, we can selectively bias the array to 


receive signals from that direction. 
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my does this selectively bias the response of the 
meeatrem tO Signals from that direction? Consider a single 


mipuese arriving from the direction of our selected beam. 
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Plotting the output of each hydrophone as a function 


Seeetime produces the following result. 


Myarophone 2 


RESPONSE 
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Hydrophone l 


RESPONSE 
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By delaying the output of hydrophone 2 for a time 
ty = At and adding the output of the two hydrophones, 
the two impulses reinforce each other and the output 


will appear as in the next drawing. 
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If, however, the impulse arrived from a different 
direction, the time difference between arrival at the two 


Mydrophones will be different. 


Hydrophone 2 


$PONnse 


Hydrophone 1 


RESPONSE 





If we input the output of hydrophone 2 through the 
seme delay, At, the plot of the summation of the two 


Signals is as shown in the next illustration. 


hie 





RESPONSE 
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Thus we see that signals arriving from the selected 
direction (time delay At) will selectively reinforce 
one another in the processor while those from other 


ferections will not. 


If we examine a 3-hydrophone linear array we observe 


that tne efiect is even more pronounced. 
ya 
i. ee 
d. z S 


By delaying the output of hydrophone 3 for a length 
@emeime 2ZAt, and the output of hydrophone Zz for a length 
of time At, a wave from the desired direction will produce 


Ie effect shown in the next figure. 
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Hydrophone 1 
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A wave from a different direction will give the 


following output. 


Hydrophone 3 


RESPONSE 


Hydrophone 2 


REAPONSE 


Hydrophone 1 


RESPOMSE 





SMeput after 
Summation 


RESTONSE 
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It 1s apparent that adding hydrophones to the array 
Mjhmes 1t more responsive to signals from the chosen 
mmmccetion. Also, since the output of the system is not 
pmeorrelation of the individual hydrophone outputs, but 
a simple summation of their energies, the system will 
pomeerOr narrow-band or single-frequency line signals 
femwel) as for broadband signals. 

By using a number of hydrophones in this manner in 
Segre tion with an energy detector, and by setting the 
@areshold high enough to preclude response to the un- 
memmeorced signals from other directions, one has a 
femyedirectional detection capability from a series of 


omnidirectional hydrophones. 
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